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ABSTRACT 


The exploration of the outer solar system using Voyagers 1 and 2 
Induced dual-wavelength radio observations during occultatlons of the 
spacecraft by Jupiter, Saturn, and Titan. One of several areas of 
research based on these observations Involves the way that refractive 
Irregularities In the atmospheres of these' planets caused amplitude and 
phase scintillations at both wavelengths (3,6 and 13 cm). 

A new technique has been developed and applied In order to 
objectively analyze the radio scintillations observed during planetary 
occultatlons of Voyager. In this method, temporal spectra of the weak 
fluctuations in anplltude and phase are Interpreted using the theory for 
scattering from a random distribution of anisotropic irregularities 
embedded In a thin screen. Comparisons between measured spectra and the 
theory produce least-squares solutions for four parameters which 
characterize the underlying atmospheric scattering structures: (1) the 

axial ratio, (2) the orientation, (3) the power-law exponent of the 
spatial spectrum, and (4) the magnitude of the spatial variations In 
refractive index. The covariance matrix of the solutions Is calculated 
by propagation of errors. 

For Jupiter and Faturn, this study considers only data pertaining 
to Irregularities In the upper ionospheres. At both planets, scattering 
was observed from Irr? ularlties of approximate size 2 to 100 km. 

Within this size interval, the three-dimensional spatial spectrum obeys 
an inverse power law; further, all observed irregularities appear to be 
anisotropic. Magnetic field orientations In the Ionospheres of Jupiter 
and Saturn were inferred from the direction of alignment of the 
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Irregularities, The occultatlon measurements conflict with the 

predictions of Jovian magnetic field models, but generally confirm 

current models of Saturn's field. 

During the occultatlon of Voyager 1 by Titan, refractive 

Irregularities of size 0.2 to 2 km distributed throughout the lowest 

90 km of the atmosphere caused appreciable signal fluctuations. 

Analysis of measured spectra yielded evidence of horizontal variations 

In gas density, and hence dynamics. Above 25 km, the altitude profile 

of Intensity scintillations closely agrees with the predictions of a 

theory based on the characteristics of freely propagating Internal 

gravity waves. These observations confirm the presence on Titan of 

atmospheric gravity waves, possibly driven by a convective region near 

the surface. The net upward flux of wave energy Is approximately 

2*10 W m -- too small to affect the thermal structure of Titan's 

lower atmosphere. An upper bound on the net upward flux of horizontal 

-3 -2 

momentum Is roughly 3*10 N m ; the transport of momentum by gravity 
waves could be Important In maintaining the zonal stratospheric 
circulation. 
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CHAPTER 1 


INTRODUCTION 


1,1 Background 

Following their launch In 1977, the two Voyager spacecraft 
successfully explored the Jovian and Saturnian systems. Kohl base and 
Penzo [1977] previously described the overall mission objectives, while 
Stone and Lane [1979a, b] and Stone and Miner [1981, 1982] summarized 
the preliminary results. Of interest here are several aspects of the 
experiments conducted through use of microwave radio signals transmitted 
by the spacecraft, and received on the earth. 

A number of properties of a planetary atmosphere can be measured by 
the radio occultation method, using the radio link between the earth and 
a spacecraft which is occulted by the planet. (For a detailed 
discussion of this measurement technique and its potential, see Eshleman 
[1973].) For this reason, the Voyager missions to the outer planets 
included radio occultations of both spacecraft by Jupiter and Saturn, 
and a Voyager 1 occultation by Titan [Eshleman et al., 1977], Through 
analysis of the dual-wavelength (3.6 and 13 cm) radio data from the 
Voyager occultations by Jupiter, Eshleman et al. [1979b, c] and Lindal 
et al. [1981] deduced altitude profiles for several quantities which 
characterize the atmosphere. These include the gas refractivity , 
molecular number density, pressure, temperature, and microwave 
absorption in the neutral atmosphere, as well as the average electron 
density in the ionosphere. Using data from subsequent Voyager 
occultations, Tyler et al. [1981a, 1982] derived analogous, though 
preliminary, results for Saturn's atmosphere and ionosphere. Analysis 
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of the effect of Titan's atmosphere on the two radio signals transmitted 
by Voyager 1 produced detailed altitude profiles for the neutral 
atmosphere, and showed that there was neither a detectable free electron 
density In the Ionosphere nor perceptible microwave absorption by 
atmospheric gases [Tyler et al.. 1981ai; Lindal et al.. 1983]. These 
results all pertain to the average structure of the planetary 
atmospheres. 

During each Voyager occultatlon, the radio receiver moved through 
the diffraction pattern associated with refractive Irregularities In the 
occulting planet's atmosphere. Rapid signal fluctuations occurred at 
both wavelengths as a result. These ubiquitous scintillations are a 
distinctive feature of the Voyager occultatlon measurements, and thus 
are Incentive to further study. Consequently, the present research was 
directed toward a primary objective -- to Infer characteristics of the 
atmospheric structures through analysis of the weakly-scattered radio 
signals. To anticipate the principal results, properties of the 
refractive irregularities were determined from the occultatlon data, and 
yielded new constraints for models of the Jovian and Saturnian magnetic 
fields, as well as a first-order characterization of small-scale 
dynamics in the atmosphere of Titan. 

Scintillations are commonly observed during occultations of 
spacecraft radio sources by planetary atmospheres [for example, see Woo 
et al., 1980a, or Woo and Yang, 1978], and are also a persistent feature 
of stellar occultations [for example, see French et al., 1983, or 
Elliot, 1979], Despite significant progress in developing a theory for 
weak scattering by atmospheric irregularities [Haugstad 1978a, b, 1979; 
Hubbard et a1., 1978], prior to this study there existed no systematic 
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and objective procedure for applying the theory toward analysis of the 
signal fluctuations. Consequently, this research was motivated by a 
second Important objective -- to establish a rigorous methodology for 
attacking this type of problem, 

1.2 Organization 

After a discussion of the refractive index of planetary atmospheres 
at microwave frequencies. Chapter 2 proposes physical models for the 
refractive irregularities In the Ionosphere and neutral atmosphere. 

These models serve to characterize the atmospheric scattering 
structures, and are used subsequently to compute theoretical 
scintillation spectra for comparison with the Voyager occultation 
measurements. In order to complete the development or the theory. 
Appendix A derives the formulas for weak scattering of microwaves from a 
slightly irregular refracting atmosphere. The results In Appendix A are 
not original, but are Included for convenience and as background for 
subsequent discussions. Appendix B establishes the limit of validity 
for the single-scattering results, and Appendix C provides a physical 
interpretation of the mathematical formulas. 

Chapter 3 di'scribes a new technique for systematical 1y analyzing 
and interpreting vhe radio scintillations observed during planetary 
occultatlons. Through use ^f conventional methods of spectral data 
analysis and the existing weak-scattering theory, this procedure yields 
least-squares solutions for parameters which characterize the 
atmospheric scattering structures. The covariance matrix of the 
solutions Is calculated by propagation of errors in order to evaluate 
the measurement uncertainties; as demonstrated in Appendix E, these 
uncertainties are realistic and not simply formal. 
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This new procedure was applied In processing the dual-wavelength 
data from the Voyager radio occultatlons by Jupiter. Saturn, and 
Titan. Chapters 4 and 5 report new results concerning the plasma 
Irregularities In the upper Ionospheres of Jupiter and Saturn. In 
Chapter 6, the measured orientation of the irregularities Is used to 
Infer the magnetic field direction at several locations In the 
Ionospheres of these two planets; the occultatlon measurements are 
compared with the predictions of current Jovian and Saturnian field 
models. On the basis of evidence derived from the Voyager 1 occultatlon 
data, Chapter 7 reports the discovery of Internal gravity waves In 
Titan's upper atmosphere. Wave parameters, Including the vertical 
fluxes of energy and momentum, are estimated, and the source of the 
gravity waves Is considered. 

Finally, Chapter 8 summarizes the results and provides suggestions 
for future research. 

1.3 A Note on Terminology 

The conventional terminology for atmospheric waves Is somewhat 
confusing; a comment on the notation used here Is therefore In order. 

In conformity with common usage In the literature, atmospheric waves in 
which gravity provides the predominant restoring force will be referred 
to as "internal gravity waves." These waves are not confined to any 
fluid Interface, but Instead propagate freely through the atmosphere. 
Gossard and Hooke [1975, pp. 13-14] discuss the terminology In detail, 
and offer the more descriptive phrase "buoyancy waves" In referring to 
this type of fluid motion. It should be emphasized that the atmospheric 
waves considered here are unrelated to the gravitational radiation 
predicted by the general theory of relativity. 
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1.4 Contributions 

The contributions by the author In the course of this research can 
be summarized as follows. 

(1) An objective procedure was developed for systematically 
analyzing the signal fluctuations observed during planetary 
occul tatlons. This new approach Is more accurate and complete 
than previous methods. 

(2) When applied to the Voyager radio occultatlon data, this 
procedure yielded new results concerning the morphology of 
Jupiter's Ionospheric Irregularities, Furthermore, the 
magnetic field orientation In the Ionosphere was Inferred from 
the alignment of anisotropic plasma structures. These 
occultatlon measurements were compared with the predictions of 
current Jovian magnetic field models; significant discrepancies 
were found. 

(3) Nety results concerning Saturn's ionosphere and magnetic field, 
analogous to those at Jupiter, were derived from the Voyager 
radio occultatlon data. However, In this case, the Inferred 
magnetic field orientations corroborate present Saturnian field 
models . 

(4) Through analysis of the Voyager 1 radio occultatlon data, 
freely-propagating internal gravity waves were discovered in 
the upper atmosphere of Titan. Gravity-wave parameters, 
including the vertical fluxes of energy and momentum, were 
estimated. These results were applied to assess the Importance 
of dynamics in determining the average thermal structure and 
circulation of Titan's atmosphere. 
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CHAPTER 2 


MODELS FOR REFRACTIVE IRREGULARITIES IN PLANETARY ATMOSPHERES 

2.1 Introduction 

During each of the Voyager radio occultatlons by Jupiter, Saturn, 
and Titan, signal fluctuations occurred at both wavelengths as a result 
of scattering from refractive Irregularities In the atmospheres of these 
planets. With the objective of constructing a theory for Interpretation 
of these data, this chapter discusses the refractive Index of 
atmospheres at microwave frequencies (Section 2.2), and proposes 
physical models for the refractive Irregularities In the Ionosphere 
(Section 2.3.1) and In the neutral atmosphere (Section 2.3.2), These 
models contain five free parameters which characterize the morphology of 
the scattering structures. When combined with the formulas In Appendix 
A, the models provide a basis for computing theoretical scintillation 
spectra. In later chapters, measured spectra are compared with the 
theory in order to determine solutions for the free parameters. 

2.2 The Refractive Index of Planetary Atmospheres at Microwave 
Frequencies 

As discussed in Appendix A, the radio frequencies used for the 
Voyager occultatlons (2295 MHz and 8415 MHz) are much larger than the 
electron gyrofrequency , the plasma frequency, and the electron collision 
frequency In the Ionosphere of Jupiter or Saturn. 


PRECEDiNO PATJE BEANK NOT FrLMF.n 

pAQfc 


7 


Models for Refractive Irregularities 


In this case, the refractive Index n(r) 1n the Ionosphere Is related to 
the electron density through 1973] 


r_ k 


n(r) = 1 - Ng(r)'Vg = 1 - Ng(r) 


( 2 . 1 ) 


Here, N„(r) Is the concentration of free electrons, V 1s the refractive 

•15 

volume of one electron, r^ = 2,82*10 m Is the classical electron 
radius, k Is the radio wavelength, and the origin for the position 
vector r Is at the center of the planet. In the expression above, a 
small correction which leads to Faraday rotation has been Ignored. As 
the Voyager spacecraft transmitted signals In right-hand circular 
polarization, this approximation Is acceptable. Note, however, that 
measurements of Faraday rotation are planned as part of an occultatlon 
experiment in a future investigation of the Jovian system using the 
Galileo orbiter, which will transmit one signal in linear polarization, 
and another in circular polarization [private communication, H. T. 
Howard, 1983]. 

In the neutral atmosphere, the refractive Index of a mixture of 
gases Including a suspension of dust particles or condensates can be 
expressed as [Eshleman, 1973] 


n(D = l*|v^> 


* Vj 

J 


( 2 . 2 ) 


where the parameters Np(r), Vp, and np are the number density, volume, 
and refractive index of the suspended particles, and the summation Is 
over the various constituents with number densities Nj(r), and constants 
of proportionality V.. For nonpolar gases, V. appears to be independent 

J J 
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of concentration, temperature, and pressure [e.g. , Tyler and Howard, 
1969], This formula remains valid under general atmospheric conditions 
when the mutual separation of the suspended particles is small compared 
to 


2.3 Physical Models for Refractive Irregularities 

The preceding section considered the principal contributions to 
refractivity in planetary ionospheres and atmospheres; here, models for 
the spatial variations in refractive index are proposed. Notice that in 
deriving theoretical expressions for the weak-scattering scintillation 
spectra in Appendix A, the refractive index in the phase screen was 
described as a locally homogeneous random field. In this case, the 
theoretical spectra (equations A36-A38 and A41-A43) depend only on the 
three-dimensional spatial spectrum of refractive index, (see equation 
A12), which constitutes a complete statistical description of the 
refractive irregularities for occultation measurements. The sections 
that follow describe models for appropriate in the ionosphere and the 
neutral atmosphere, and define the variables and coordinate systems that 
are used in subsequent chapters. For a discussion of the physical 
significance of 4^, the reader is referred to Tartarskii [1971, 
pp. 1-46] or Yaglom [1962]. 

2,3,1 Ionospheric Irregularities. It is useful to consider the 
terrestrial ionosphere in order to construct a theoretical model for the 
ionospheric irregularities at Jupiter and Saturn. In the earth's upper 
ionosphere, the spatial spectrum of irregularities in electron density 
obeys a power law for length scales ranging from the ion gyroradius (a 
few meters) to the ionospheric scale height (a few tens of kilometers) 
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[Yeh and Liu, 1982; Booker, 1979; Rufenach, 1972, 1974; Dyson et al,, 
1974; Phelps and Sagalyn, 1976], In addition, the Irregularities are 
anisotropic and "field aligned" -- the correlation length of these 
structures Is largest along the direction of the earth's magnetic field 
[Spencer, 1955; Singleton, 1970a; MacOougall, 1981], By analogy, the 
Irregularities at Jupiter and Saturn are expected to possess these same 
characteristics for two reasons. First, the magnetic field In the 
terrestrial Ionosphere [Langel et al,, 1980] Is comparable In magnitude 
to Saturn's field [Connerney et al., 1982], but Is greatly exceeded by 
Jupiter's field [Acuna and Ness, 1976]. Hence, the magnetic field is 
expected to constrain the motion of charged particles In the uppe*" 
ionospheres of these two planets In the same manner as It does at the 
earth. Second, the size of the Irregularities observed at Jupiter and 
Saturn (see Chapters 4 and 5) was well within the bounds of the Ion 
gyroradlus (* 1 meter) and the Ionospheric scale height (« 10 km) 
[Eshleman et al., 1979b, c; Tyler et al., 1981a, 1982], Consequently, 
an anisotropic (field-aligned), power-law spatial spectrum was used in 
computing theoretical scintillation spectra for comparison with the 
measurements at Jupiter and Saturn. In a Cartesian coordinate system 
aligned with the principal axes of the irregularities, this spatial 
spectrum is expressed as [cf. Singleton, 1970b; Rino and Fremouw, 
1977] 



where Q is a normalization constant that Is proportional to the vfjriance 
of refractive Index, <n^^>; 2 h/x:q corresponds to the outer scale; 
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<2* are the components of the wave vector: d, e, and f are 

numerical factors that determine the degree of anisotropy; p is the 
power-law exponent; and the ic2-direction is parallel to the magnetic 
field, L 

For simplicity, it is assumed that the irregularities are 

symetrical with respect to rotations about the magnetic field 

direction, so that d => f. In this case, the correlation lengths of the 

ionospheric irregularities along the magnetic field and perpendicular to 

this direction are in the ratio e/f [cf. Singleton, 1970b], Large axial 

ratios (e/f » 1) correspond to highly elongated field-aligned 

structures, smaller axial ratios characterize irregularities with 

ellipsoidal shapes, and e/f = 1 for isotropic irregularities. With the 

proper choice of the principa! axis, the i rregularities can always be 

represented by an axial ratio e/f > 1. 

The appropriate scattering theory for the ionospheric occultations 

considered here is the one derived in Sections A. 2 and A. 3, where the 

average refractive index in the phase screen is homogeneous; for present 

purposes differential refraction had a negligible effect. In computing 

the theoretical spectra (equations A36-A38), must be expressed m the 

(k , K , K ) coordinate system defined in Section A.4 and shown in 
X y z 

Figure 2.1. In this system, < is parallel to the direction of radio- 

X 

wave propagation, and tc is parallel to the projection of the magnetic 

•J 

field onto the phase screen. The transformation gives 


S (k =0, K , K ) = 

n - X ’ y’ z' 


[1 + f 


(2.4) 


As shown in Appendix A, refractive variations along the direction of 
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propagation do not affect the signal scattered at small angles; 
hence, is set equal to zero in the argument of (see equations ASS- 
ASS). Rino and Fremouw [1977] previously used a power spectrum of this 
form to describe irregularities In the earth's ionosphere. 

The parameter e' is related to e, f, and the angle, A, between 
the B and the phase screen through e' = e cos (A) + f sin (A) , 

However, the effect of A is neglected here, since for the Voyager 
occultati?ns by Jupiter and Saturn it is so small that e* “ e. 

For a phase screen with a homogeneous average refractive index, 
only the relative orientations of the raypath velocity in the phase 
screen, v*, and the principal axis of the irregularities affect the 
scattering (cf. Appendix A). Hence, the phase screen coordinates 
defined in Figure 2.1 are convenient for representing the 
irregularities: a single angle, y, specifies the unknown orientation of 

the principal axis, < , relative to v*. A total of five free parameters 

•J 

(Q, p, e/f, Kq, and y) constitute the theoretical description of the 
ionospheric irregularities. 
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Figure 2.1 Schematic of phase screen for the case cf scattering 

from ionospheric irregularities (Sections A. 2 and. A.3), In this 

coordinate system, is parallel to the direction of radio-wave 

propagation and perpendicular to the screen. The magnetic field, 5, 

intersects the screen at angle A, and is parallel to the projection 

of B onto the screen. The angle y specifies the orientation of the 

principal axis of the irregularities, < , relative to the raypath 
^ y 

velocity, note that v* is perpendicular to < , From the point of 
view of a terrestrial observer, v*, k , and < are in the plane of the 
sky. 
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2.3.2 Refractive Irregularities In the Neutral Atmosphere, In 
order to Interpret the scintillations observed during the occultatlon of 
Voyager 1 by Titan, the spectrum of refractive Irregularities in the 
neutral atmosphere is modeled using a formula analogous to equation 
2.3. The motivation for choosing a power-law spectrum comes from 
observations of the earth’s atmosphere [Vinnichenki , 1970] and 
theoretical considerations [Tatarskii, 1971, p. 46-102]. This type of 
spectrum Is believed to be valid under many circumstances for describing 
the Irregularities that result from either turbulence or a superposition 
of Internal gravity waves [Dewan, 1979]. 

To account for the Important Influence of Titan's finite 
atmospheric scale height, the scattering theory must include' the effect 
of Fresnel zone distortion; consequently, the formulas derived In 
Sections A. 5 and A. 6 (equations A41-A43) are appropriate for these 
calculations. Some slight modification of from the previous form Is 
needed. In this case, must be expressed In the k^) 

coordinates defined In Section A. 4 and shown 1n Figure 2,2, where is 

parallel to the direction of radio-wave propagation, but is now 
parallel to the visible limb of the planet at the occultatlon point: 


<6 (k = 0, 

n^ X * 




11 * 


-2 


(Ak^ 


Bic K + 

y z 


Ck/)J 




. (2.5) 


The parameters A, B, and C are related to e/f and the angle 4 between 
the principal axis of the Irregularities and through 
A = e^crjs”(<P) + f^s1n^(<P) , B = 2 (e^- f^) sin(4») cos(4>) , and 
C = e sin (<i>) + f cos (<!') . As mentioned above, these occultatlon 
measurements are Insensitive to refractive variations along the 
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direction of propagation. Consequently, a terrestrial observer can 
discern only tne tl:? orientation of the irregularities in the plane of 
the sky {y-z plane) as specified by 



Figure 2.2 Schematic of phase screen for scattering from 
i rregul arities in the neutral atmosphere {cf. Sections A. 5 and A. 6). In 
this coordinate system, is parallel to the direction of propagation 
and perpendicular to the screen; ic^ is parallel to the visible limb of 
the planet at the occultatlon point; and k coincides with the vertical 
direction, opposite to the gravitational acceleration, g. The 
angle y specifies the known orientation of the raypath velocity, v*, 
with respect to the limb of the planet, Note that v* and g are 
perpendicular to 
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Figure 2.3 shows a diagram of "typical" refractive Irregularities 
as viewed from earth. The two characteristic dimensions of the 
structures (parallel and perpendicular to the direction Indicated by <|») 
are In the ratio e/f. For example, when < 1 ^ = 0, large axial ratios 
(e/f >> 1 ) correspond to layer-like structures with little horizontal 
Violation, smaller axial ratios (e/f > 1) characterize Irregularities 
with both horizontal and vertical variations, and e/f = 1 for Isotropic 
Irregularities. 



Figure 2.3 Typical irregular structures in a planetary atmosphere 
from the point of view of a terrestrial observer. The angle 41 gives the 
orientation of the principal axis with respect to the visible limb of 
the planet. 
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As 1n the representation of Ionospheric Irregularities, the free 
parameters in the model for In the neutral atmosphere include Q, p, 
e/f, and However, for scattering from the neutral atmosphere, where 
Fresnel zone distortion Is appreciable, both the relative and absolute 
orientations of v* and the principal axis affect the scattering, so that 
two angles are required to specify the geometry. In this case, y 
indicates the known direction of v*, while <i» specifies the unknown 
orientation of the principal axis; both angles are measured relative 
to Ky. These five free parameters (Q, p, e/f, tc^, and 4< ) constitute 
the theoretical description of Irregularities In the neutral atmosphere. 
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CHAPTER 3 


DATA ANALYSIS AND INTERPRETATION 


3.1 Introduction 

This chapter describes a new procedure for the systematic analysis 
and Interpretation of radio scintillations observed during planetary 
occultatlons. Through the use of conventional methods for spectral data 
analysis along with the wel 1 -developed weak-scattering theory, this 
procedure yields results which are more accurate and complete than those 
reported previously (see Section 4.1 for a comparison with the method 
used to analyze data from the occultatlons of the Pioneer spacecraft by 
Jupiter). The present technique, which employs a rigorous scattering 
theory, has evolved from the approach developed by Tyler et al. [1981b], 
where an approximate theory was applied to Interpret observations of 
radio-wave scattering from the solar corona. 

3.2 Fundamental Measurements 

The radio signals used for the Voyager occultatlon measurements 
were derived from a highly stable, radiation shielded, quartz oscillator 
carried by the spacecraft. Two unmodulated signals were transmitted in 
right-hand circular polarization through the 3.7-m spacecraft antenna, 
one at 13-cm wavelength (26 W) and a second at 3.6-cm wavelength (12 W), 
and received with the NASA Deep Space Network (DSN) 64-m antenna 
tracking stations. The spacecraft and other radio equipment used In 
this experiment has been described in more detail elsewhere [Eshleman et 
al., 1977]. 
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Analysis and Interpretation 


The DSN tracking st<^tions coherently sampled and recorded the 
signals received from Voyager for later processing. After delivery to 
Stanford, the data were first digitally filtered to bandwidths of about 
I kHz. Then, a previously developed method [Lipa and Tyler, 1979] was 
applied to estimate the Instantaneous frequency, f(t), and amplitude, 
A(t), of each received signal from a sequence of power spectra of the 
filtered data. The measurements of A(t) are unbiased in that the 
procedure accounts for the average effect of system noise. Values of 
f(t) and A(t) were typically computed at intervals. 6t = 0.05 s, although 
the sampling rates varied among occultations as a result of differences 
in the recording bandwidth of the DSN equipment and the occultation 
geometry. The specific sampling rates will be given in connection with 
the discussions in Chapters 4, 5, and 7. All subsequent processing 
steps were based on these measured values of f(t) and A(t). 

To illustrate certain signal characteristics, it is useful to 
consider the data concerning Jupiter's ionosphere which were obtained at 
Voyager 1 occultation immersion. The time history of scintillations can 
be displayed on an Argand diagram by plotting the normalized phasors 
[A(t)/<A>]“exp[j4>(t)3 in the complex plane. The normalizing value <A> 
represents an average over time, while the fluctuations in signal 
phase i>(t) were obtained by fitting a quadratic curve to measurements of 
f(t) over periods of 6.6 s and integrating the residuals. The quadratic 
detrending of f(t) before integration removes the effects of spacecraft 
motion relative to the DSN station from il>(t). In the absence of noise, 
a signal that had propagated through an ionosphere without refractive 
irregularities would appear as a constant unit vector along the positive 
real axis. 
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Figure 3.1A shows such a plot of 128 measurements of amplitude and 
phase, at 0.051 s Intervals, for data preceding Voyager 1 occultatlon 
Immersion. The fluctuations are Indicative of the noise level of the 
occultatlon measurements. At the sampl Ing rates considered here,. 
thermal efects at the DSN station and In the earth's atmosphere were the 
predominant source of amplitude and phase noise at 13 cm, and amplitude 
noise at 3.6 cm. In addition to the thermal noise sources. Instrumental 
effects — possibly In the DSN programmed oscillator or the spacecraft 
oscillator -- contributed appreciably to the phase noise at 3.6 cm. 

Figure 3, IB shows a second plot for data at a later time as the 
path of propagation from spacecraft to earth traversed Jupiter's upper 
ionosphere. The Increased but still small fluctuations In amplitude and 
phase at both wavelengths are the result of scattering from refractive 
Irregularities. The effect at the 13-cm wavelength Is larger, primarily 
as a result of plasma dispersion. Figure 3. 1C shows r third plot 
obtained at a slightly later time, and hence with ray perlapsis deeper 
in Jupiter's Ionosphere, where stronger spatial variations In free 
electron density produced a pronounced effect at 13 cm as indicated by 
the random walk of the (A, 4>} phasor about the origin. The 3.6-cm 
signal remained weakly perturbed during the corresponding time. 
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Analysis and Interpretation 


Appendix B sets a limit for the validity of the weak-scattering 
theory. When this criterion Is applied to the results In Figure 3.1, 
the data at both wavelengths in part B, and at 3.6 cm In part C, are 
found to satisfy the condition for weak scattering, while the 13-cm data 
In part C exceed this limit and correspond to strong scattering. In 
general, for weakly-scattered signals, the phasors remain within a 
relatively small region In the complex plane;; for strong scattering, the 
distribution of signal phasors has a "random" appearance. In this 
regard, Argand diagrams provide a useful diagnostic tool for evaluating 
the signal characteristics. The reader is referred to Flatte [1979] for 
additional discussion of the Interpretation of phasor diagrams, and to 
Hinson and Tyler [1983] and Hinson [1983] for comparison with 
occultation results from Titan's neutral atmosphere and Saturn's 
Ionosphere, respectively . 


3.3 Spectral Data Analysis 

The scintillation data were spectrally analyzed over time Intervals 
In which the fluctuations In amplitude and phase appeared to be 
statistically similar, or stationary, as determined by Inspection of 
Argand diagrams and the time history of A{t) and f(t). The fluctuations 
in amplitude were normalized according to 


<A> 


(3.1) 


and values of <l)(t) were computed as before using the same detrending 
interval as for determining <A>. 



Analysis and Interpretation 


Discrete sample spectra were computed from the formulas [cf, Jenkins and 
Watts, 1968, pp. 209-211, 340-344], 






I a(t„) ] 

n«0 


2 


C 



m-1 

i ♦(t„) 
n«0 


2 


(3.2a) 


and using the a(t) and 4>(t) measured at the same radio wavelength, 

, m-1 m-1 

(3.2b) 

1 / 

where f|^° ^5^73^ , k >= 0, 1, 2, ... , m-1 


Here, f. and t represent the m sampling frequencies and times in the 
K n 

discrete Fourier transform pair, and vertical bars denote the 
magnitude. In order to reduce the bias from detrending, the detrending 
interval was selected to be 2m'*6t -- twice as long as the time interval 
for each spectrum. Several sample spectra from within a time interval 
were averaged to produce smoothed spectra for amplitude, ^a^(^l^)» snd 
phase, ^^^(^l^)* and the smoothed cross amplitude-phase spectrum, 

C’g^(fk). The time interval chosen for a sample spectrum resulted from a 
compromise between long time intervals, which reduce bias effects and 
improve frequency resolution, and short intervals, which allow more 
smoothing and reduce the variance [Jenkins and Watts, 1968, pp. 230- 
257]. Typically, eight sample spectra, each computed from m = 128 time 
domain samples, were averaged to produce each smoothed result. 
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The cross spectrum provides a measure of the relative size and 
phase of correlated fluctuations at each frequency In two stochastic 
processes. Of particular Interest is the smoothed squared coherency 
spectrum defined as 








(3.3) 


The squared coherency plays the role of a correlation coefficient 
defined at each sample frequency fj^, and Is used below along with the 
phase of to demonstrate the consistency of the data with weak- 
scatterlng theory. 

analyzing the scintillations observed during the occultatlon of 
Voyager 1 by Titan's neutral atmosphere, two slight deviations from the 
above procedures were Introduced. First, as It Is conventional to 
discuss the Intensity, I(t), rather than A(t) during occultatlons by 
neutral atmospheres, the spectrum of Intensity, Instead of was 
measured. This change Is Insignificant for the case of weak scattering, 
where Tjj = 4*^ , (see Section A. 7). Second, the normalized Intensity 
was computed from 


1(t) 


. I(t) - <I(t)> 


(3.4) 


Here, the normalization <I(t)> represents a least-squares quadratic fit 
to I(t) over the time period used, and accounts for the gradual changes 
in I(t) resulting from differential refraction (see Appendix D) and 
antenna pointing, which were insignificant during the ionospheric 

2 

occultatlons. In addition, the fractional variance of Intensity, , 
was calculated for data from the Titan occultatlon. 
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3.4 Interpretation of Scintillation Spectra 

Smoothed amplitude and phase spectra, obtained by the procedure 
described In Section 3.3, were used as estimates of the actual spectrum 
of the data. To a good approximation, the measurements at the discrete 
sample frequencies are uncorrelated and obey chi-squared statistics 
[Jenkins and Watts, 1968, pp. 230-2573, with typically 15 degrees of 
freedom for the data analyzed here. Noting that the chi-square 
approaches a Gaussian distribution with Increasing degrees of freedom, 
the statistics of the measurements can conveniently and accurately be 
modeled as Gaussian with fractional variance where M is the number 
of sample spectra averaged In the smoothed result. 

In order to identify the contribution to signal fluctuations that 

arises from atmospheric scattering, the noise level of the spectra was 

determined empirically from measurements before Immersion and after 

emersion, Thei.e noise levels are Indicated, for example, 1n Figures 

4.3, 5.5, and 7.3, below; their appearance can be explained as 

follows. In the absence of atmosplierlc scattering, successive temporal 

measurements of a{t) or f(t) are uncorrelated; hence, the spectra of 

noise fluctuations 1n both a(t) and f(t) are white. By applying basic 

Fourier transform properties to this result, the "power spectrum" of 

noise fluctuations In <{>(t} has slope -2 on a logarithmic display, 

2 

All spectra analyzed in Chapters 4, 5, and 7 exhibit <a > < 0.06, 
which satisfies the limit for validity of the weak-scattering theory 
(Appendix B). In order to interpret the weak ionospheric 
scintillations, theoretical spectra were computed numerically from 
equations A36-A38 using the expression for in equation 2.4 as a model 
for the refractive Irregularities. Alternatively , theoretical spectra 
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for weak scattering from the neutral atmosphere were computed 
numerically from equations A41>A43 and 2.5. In the models 
for (equations 2.4 and 2.5 In Section 2.3), the free parameters which 
characterize the refractive i rregularltles Include: (1) the Inverse 

power-law exponent, p; (2) the axial ratio, e/f; (3) the orientation, 
specified by y In the Ionosphere, or <1^ In the neutral atmosphere; and 
(4) the strength, Q. (The fifth parameter, k^. Is discussed below.) 

The theoretical spectra were fitted to each measured amplitude or phase 
spectrum In a least-squares sense by varying these four free parameters 
to minimize the residual 

l»9ioLf(fk)J - • (3.5) 

Here, ^(tj^) is the measured spectrum, w(f|^) Is the theoretical spectrum, 
and the summation runs over sample frequencies, f|^, for which the 
scintillation spectra exceeded the noise level by at least 3 dB. In 
some cases computational convenience led to larger values for this 
cutoff. Notice that the spectra are compared in the log domain, where 
the measurements at all sample frequencies have the same variance. The 
values of the parameters at the minimum were taken as characteristic of 
the irregularities, which requires a homogeneous over the atmospheric 
region corresponding to one spectrum. 

As a test of these solutions for parameters, spectra were analyzed 
by first subtracting the noise level and then fitting the theoretical 
spectra. The parameter estimates were found to be insensitive to this 
change in procedures, and hence to the measurement noise. 
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Within the limitations Imposed by experimental geometry and the 
method of data analysis, no evidence for an outer scale, 2w/tc^, was 
found In any case. Consequently, all theoretical spectra were 
calculated from the formulas of Chapter 2 and Appendix A In the limit 
where 0; this parameter then had a negligible effect on the theory 
over the range of measurement frequencies. 

The variance of refractive Index was determined from (see equation 

A12) 


= /]/ *n('^x’^y’''z^ ^*^x ^“^y 

where n^(r) represents the variations In refractive Index away from the 

spatial average, and the parameters of were determined from the 

data. The limits of Integration were restricted to the range of wave 

vectors which could be detected, so measured values of <n^^> given below 

apply only to this range. Using the formulas In Section 2,2 (equations 
2 

2.1 and 2.2), <n^ > can be related to variations In free electron 
concentration In the ionosphere or fluctuations In particle densities in 
the neutral atmosphere, as appropriate. 

The covariance matrix of the solutions for the four parameters, 
[Cov], was computed by standard propagation of errors [Brandt, 1970, 
pp. 196-212] using the theoretical expressions in t;quat1ons A36 and A37 
or A41 and A42 as models for the process. The general form of the 
matrix Is 

[Cov] = P]"^ . ( 3 . 7 ) 
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Here, P Is the tn-by-4 matrix of first partial derivatives of the 
theoretical spectrum with respect lo the four parameters at each of the 
m sample frequencies. The variance of each measured parameter Is 
proportional to {see equation 3.5). As discussed by Brandt, the 
diagonal elements of [Cov] can be considered the variance of the four 
parameter measurements. While the off-diagonal elements are useful In 
specifying the statistical relationships among the parameters, they are 
not reported here as they do not give additional Information on the 
uncertainty In any one parameter. All uncertainties given In tables In 
Chapters 4-7 refer to the standard deviations as calculated from 
equation 3.7. On the basis of the discussion In Appendix E, these 
uncertainties are believed to be realistic and not simply formal. 

As a test of this procedure for computing uncertainties, the 
measurements were simulated using a statistical model. Artificial 
spectra were generated by adding different realizations of uncorrelated 
Gaussian noise with fixed fractional variance to a theoretical 
spectrum. Then, each of these spectra was analyzed by the procedure 
described above to obtain least-squares solutions for parameters. The 
statistics of the parameter measurements obtained in this way confirmed 
the predictions of the covariance matrix. 

The smoothed cross spectrum, provides an estimate of the true 
cross spectrum of two stochastic processes. The variance of the 
estimate can be reduced through averaging, but also depends on the 
uncontrollable influence of the coherency spectrum. At frequencies 
where the coherency Is small, the processes are only weakly correlated, 
and the estimates of are inherently noisy. This effect resulted In 
a large variance in for the present measurements of weak 
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scattering. Hence, the solutions for parameters obtained from analysis 
of were less reliable than solutions derived from or and 
will not be given here. However, for completeness and as a test of the 
theory. Chapters 4 and 5 Include comparisons of theoretical spectra with 
the experimentally derived squared coherency and phase of *'a<i)* 
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CHAPTER 4 


SPATIAL IRREGULARITIES IN JUPITER'S UPPER IONOSPHERE 
4,1 Introduction 

Following the firat two spacecraft encounters with the Jovian 
system by Pioneer 10 In December, 1973, and Pioneer 11 In December, 1974 
[Hall, 1975; Gehrels, 1976], Voyagers 1 and 2 were used to continue the 
reconnaissance of Jupiter and Its environment during flybys In March and 
July, 1979 [Stone and Lane, 1979a, b]. The Investigations with each of 
the four spacecraft Included radio occultatlon measurements of Jupiter's 
Ionosphere and neutral atmosphere. Through analysis of the radio data 
from the occultatlons of Voyagers 1 and 2, Eshleman et al. [1979b, c] 
and Linda! et al. [1981] determined vertical profiles of the gas 
refractlvlty , molecular number density, pressure, temperature, and 
microwave absorption In the neutral atmosphere, along with the average 
electron density In the Ionosphere. Results of a similar nature, as 
derived from the Pioneer radio occultatlon data, were reported by Kliore 
and Wolceshyn [1976] and Fjeldbo et al. [1976], In the present 
discussion, attention Is focused on rapid fluctuations of the radio 
signals observed during the Voyager occultatlons by Jupiter's 
1 onosphere. 

During the occultatlons of Pioneers 10 and 11 by Jupiter, 
ionospheric structures caused scintillations in the radio signal 
transmitted by the spacecraft and received on the earth. Woo and 
Yang [1976, 1978] and Woo and Armstrong [1980] analyzed these data and 
attributed most of the amplitude fluctuations to scattering by 
anisotropic turbulence. The spectrum of the turbulence appeared to 
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follow a power law with exponent near the Kolmogorov value of -11/3 [Woo 
and Yang, 1976, 1978], although subsequent analysis of spectral 
broadening suggested a value of -3.3 [Woo and Armstrong, 1980]. In 
addition, the magnetic field orientation In Jupiter's Ionosphere was 
Inferred from the measured direction of anisotropy. At one location, 
the estimate of field orientation agreed with the predictions of a 
Jovian magnetic field model derived from Pioneer magnetometer data; at 
three other locations the measurements disagreed with the model by as 
much as 10® [Woo and Yang, 1978; Smith et a1., 1976], However, Woo and 
Yang did not Indicate the uncertainties of the measurements, so the 
significance of the differences Is unclear. 

When comparing the Pioneer and Voyager results concerning 
Ionospheric Irregularities, It should be noted that the Voyager data are 
Improved with respect to the Pioneer data for three reasons. First, 
Voyager observations were obtained simultaneously at two coherently 
related radio wavelengths (3,6 and 13 cm) instead of the single 13-cm 
wavelength of the Pioneer data. For the same distribution of plasma 
Irregularities, the 3.6-cm signal Is much less susceptible to scattering 
than the 13-cm signal -- the variance of plasma induced intensity 
fluctuations was typically 30 times larger at the longer wavelength. 
Consequently, the 3.6-cm signal was weakly perturbed In regions of the 
lower Ionosphere where the 13-cm signal experienced strong scattering. 
The short-wavelength signal thereby significantly extended the range of 
the occultatlon measurements, since weak scattering Is well understood, 
whereas strong scattering Is difficult to Interpret, Second, for 
Voyager at 3.6 cm the signal-to-nolse ratio (SNR) was about 100 times 
greater than for Pioneers 10 and 11. The Improved signal revealed more 
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detail In the measured spectra, which allowed closer comparison between 
theory and observations (compare Figures 7-10 of Woo and Yang [1978] 
with Figures 4,3, 4.6, and 4.10 In this chapter). Third, the Increased 
SNR combined with the Improved stability of the Voyager radio equipment 
allowed analysis of fluctuations in signal phase, which was not possible 
with Pioneer. 

There also has been a significant advance In the methods used for 
analyzing scintillation spectra (Chapter 3). This study established a 
consistent procedur-e for evaluating the quality and character of the 
occultatlon signals (e.g., Argand diagrams, Section 3.2). Furthermore, 
four parameters of the theory were determined by the method of least 
squares, and the measurement uncertainties were calculated by 
propagation of errors (Section 3.4). Woo and Yang [1976, 1978] did not 
report the use of an objective method in obtaining a fit of the theory 
to the measured Pioneer spectra. In addition, they frequently showed 
the same theoretical parameters In fits of theory to spectra measured at 
different altitudes (see Woo and Yang [1'378], Figures 7-10), which 
leaves some question as to whether the best fit has been found. 

Finally, as the uncertainties in final parameters were not given, the 
Pioneer results were dif'icult to Interpret. 

Section 4.2 describes the Voyager occultatlon measurements; Section 
4.3 follows with a discussion of findings for the spatial 
characteristics of the ionospheric irregularities. Comparisons with the 
previously published Pioneer results are included where appropriate. 
Present results, which have been published elsewhere [Hinson and Tyler, 
1982], are used subsequently 1n Chapter 6 to infer magnetic field 
orientations in Jupiter's ionosphere. 
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4.2 Experimental Geometry and Measurements 

Figure 4.1 shows the view from earth of the Voyager occultations 
including the two raypath trajectories in the atmosphere. For Voyager 1 
(VI), occultation immersion occurred on March 5, 1979, at 10®S, 63®W, 
with solar zenith angle 6 » 82“ (late afternoon). Occultation emersion 
followed at 0“, 314“W, 6 = 98“ (predawn). Voyager Z (V2) occultation 
immersion occurred on July 10, 1979, at 65“S, 260“W, 6 >= 88“ (evening), 
with emersion at 50“S; 150“W, 6 = 92“ (early morning). All longitudes 
above are expressed in System III, 1965 [Seidelmann and Divine, 1977], 
Prior to immersion and immediately following emersion, the radio path 
from the spacecraft to earth traversed only Jupiter's ionosphere; at 
these times scintillations of both radio signals were observed at all 
four locations. 


N 



Figure 4.1 Schematic view from earth of the Voyager occultations 
by Jupiter showing the paths followed by the spacecraft radio Images- 
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Several parameters which characterize the Voyager occultations by 
Jupiter near immersion and emersion are listed in Table 4.1. Here, 

V* is the raypath velocity in the plane of the sky, A specifies the 
orientation of v* with respect to Jupiter (see Figure 6.1), D is the 
distance from the spacecraft to the visible limb of Jupiter, and the 
Fresnel zone radius is given by where K denotes the radio 

wavelength. 


Table 4.1 Geometric parameters that characterize the Voyager 
occultations by Jupiter near immersion and emersion. 


Space= 

craft 

Lat.^ 

p 

LonQo^ 

!vj» 

km/s 

A, 

degrees 

Fresnel-zone radius, km 
0. km 13 cm 3.6 cm 

VI 

10“S 

63“W 

20.1 

175. 

4.45«10^ 

7.62 

3.98 

VI 

0” 

314“W 

20.2 

5. 

5.47-10® 

8.46 

4.42 

V2 

65®S 

260"W 

10.5 

122. 

1.35-10® 

13.3 

6.94 

V2 

50“S 

150<»W 

10.5 

57. 

1.42-10® 

13.6 

7.13 


^lanetocentric 

^Sys III, 1965 [Seidelmann and Divine, 1977] 


The procedure for preliminary data reduction was described in 
Section 3.2. At both wavelengths, values of signal frequency, f(t), and 
amplitude, A{t), were determined each 0.051 s for VI and 0.102 s for V2 
in order that signal parameters would be adequately sampled. 

Limitations of the radio equipment and observation time resulted in 
uncertainties in a single estimate of frequency and amplitude of 
approximuLely 0.2 Hertz and 2%, respectively. The processing steps in 
the remainder of this chapter, as well as the Argand diagrams in Figure 
3.1, were based on these measured values of f(t) and A(t). 


35 


Jupiter's Ionosphere 


Using the method described In Section 3.3, the fluctuations In 
phase <ti(t) and normalized amplitude a(t) were reduced to obtain temporal 
spectra. Theoretical spectra were then fitted to the measurements by 
the least-squares procedure In Section 3.4, yielding solutions for the 
four principal parameters that characterize the ionospheric 
irregularities (cf. Section 2.3.1). The parameter measurements and 
their uncertainties are discussed below. 
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4.3 Results 

VI Immersion (IQ”S, 63”M). The profile of the average electron 
density In Jupiter's upper Ionosphere appears In Figure 4.2. This curve 
was derived from the measurements during VI Immersion [Eshleman et a1., 
1979b], and can be extended to lower altitudes using more elaborate 
inversion techniques [Fjeldbo, 1964]. 



Figure 4.2 Profile of average electron concentration In the upper 
Ionosphere at VI Immersion. Regions 1-4 indicate altitude intervals for 
computing spectra. These altitudes apply to periapsis of the radio 
raypath from spacecraft to earth, and are referenced to the 1 mbar level 
in the neutral atmosphere [profile from Eshleman et al., 1979b]. 


Jupiter's Ionosphere 


Weak scattering of the 3.6-cm signal was observed throughout the 
Ionospheric occultation; spectra of 3.6-cin amplitude and phase from the 
altitude Intervals 2, 3, and 4 of Figure 4.2 appear In Figure 4.3. The 
13-cm signal experienced strong scattering 1n the lower Ionosphere; 
since analysis here Is restricted to weak scattering, a 13-cm spectrum 
was computed only over a single high altitude Interval, region 1 of 
Figure 4.2. In the spectra, the Indicated noise levels were derived 
from measurements before Immersion (see Section 3.4). No phase spectra 
are shown for the high altitude smoothing intervals, regions 1 and 2, 
because noise dominated these measurements. The best fit theoretical 
spectra computed from equations A36, A37, and 2.4 are shown superimposed 
In Figure 4.'^ . 

Figure 4.4 displays the observed squared coherency and cross 
spectrum phase for the two smoothing Intervals where both the phase and 
amplitude spectra were available (regions 3 and 4), The 1-a 
uncertainties depend on the squared coherency and are therefore a 
function of frequency — one typical value Is indicated. The high 
frequency portion of each spectrum corresponds to measurement noise for 
which the cross spectrum phase has a random distribution. Using the 
parameters obtained from analysis of the amplitude and phase spectra 
(Figure 4.3), the theoretical squared coherency and cross spectrum phase 
were computed and superimposed on each measured spectrum. The agreement 
between theory and observations -- especially at the striking shifts In 
the phase of the cross spectra — supports the validity of the parameter 
measurements and the underlying theory. 
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PBEQUENCY. Hi FREO'.IENCV, Hi 

Figure 4.3 Amplitude and phase spectra at 3.6 and 13 cm from the 
smoothing intervals (regions 1-4) In Figure 4,2. The I -0 uncertainties 
in each panel apply at all sample frequencies. Theoretical spectra 
(smooth curves superimposed) were fitted In a least-squares sense to 
measured spectra over the frequency range where the amplitude or phase 
fluctuations due to Ionospheric scattering exceeded the noise levels by 
at least 3 dB. 
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FREQUENCY, Hi FREQUENCY, Hz 


Figure 4.4 Squared coherency and cross spectrun phase at 3.6 cm 
from regions 3 and 4 of Figure 4.2. The l-o uncertainties are indicated 
at one sample frequency. The high frequency portion of each spectrum 
corresponds to measurement noise for which the cross spectrum phase has 
a random distribution. Theoretical curves are superimposed. 
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The least-squares solutions for Ionospheric parameters derived from 
amplitude and phase spectra are listed In Table 4.2, with the l-o 
uncertainties as calculated by propagation of errors (cf. Section 
3.4). The altitude assigned to each measurement corresponds to the 
center of the appropriate smoothing Interval; each measurement should be 
Interpreted as a spatial average over this altitude range. As indicated 
by the close agreement between theory and measurements in Figure 4.3, 
the three-dimensional spatial spectrum. closely follows an Inverse 
power law with exponent, p, between 3.3 and 3.7. It Is difficult to 
discern any vertical variation In p because of the measurement 
uncertainties; however, the axial ratio, e/f, apparently varies with 
altitude. Strong anisotropy was detected high In the topside ionosphere 
where collision frequencies are low -- the magnetic field Is thought to 
strongly constrain the motion of charged particles in this region. 
Randomizing collisions occur more frequently in the lower ionosphere and 
reduce magnetic field effects; this may explain the decrease In the 
axial ratios at lower altitudes. The six estimates of y, the 
orientation of Jupiter's magnetic field relative to the raypath velocity 
(see Figure 2.1), are mutually consistent and show no variation with 
altitude. This was expected, since models of Jupiter's field predict no 
appreciable variation with radius over the extent of the upper 
ionosphere. 

In the weak-scattering theory applied here. Ionospheric 
Irregularities are assumed to obey the Taylor hypothesis so that they 
can be considered "frozen" on time scales of interest. The only 
significant time variation results from motion of the radio path across 
a fixed field of Irregularities. For highly anisotropic irregularities. 
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only the component of ray path motion of magnitude |v^|»s1n(Y) which Is 
perpendicular to the elongated axis will result In signal scintillations 
{cf« Appendix C). On the basis of this Interpretation, Irregularities 
of approximate size 2 to 100 km contribute to scattering over the 
frequency interval 0.2 to 10 Hz In Figure 4.3. These lower and upper 
limits result froin noise limitations and the length of the measurement 
time Interval, respectively, and do not represent either Inner or outer 
scales. The spatial fluctuat/ons of elecron density, were 
determined from the parameter Q (Section 2.3.1) by using equatloris 3.6 
and 2.1. The values listed In Table 4.2 refer )nly to Irregularities In 
the stated size Interval, and correspond to fractional variations In 
free electron concentration of about 10 to 15%. 


Table 4.2 Characteristics of Irregularities In Jupiter's upper 
Ionosphere at VI Immersion. 


Height, Spectrum Y, Power-law Axial 6N , cm 

km type degrees index ratio 


1200 

3.6-cm 

a(t) 

56 

+ 

2.4 

3.6 

± 

0.11 

4.5 


0,3 

14000 

+ 

3000 

1200 

3,6-cm 

f(t) 

56 

± 

1.8 

3.4 

± 

0.06 

4.5 

± 

0.8 

16000 

± 

3000 

2200 

3.6-cm 

a(t) 

SB 

+ 

1.7 

3.7 

± 

0.10 

6.0 

± 

1.3 

12000 

t 

2000 

2200 

3,6-cm 

1>(t) 

58 

± 

1.7 

3.5 

± 

0.06 

4.5 

± 

0.8 

12000 

i 

2000 

3200 

3,6-cm 

a(t) 

56 

± 

1.0 

3.6 

± 

0.12 

10. 

± 

2.2 

2700 

t 

500 

3600 

13 -cm 

a(t) 

58 

i 

1.6 

3,3 

± 

0.08 

7.5 

± 

1,7 

1400 

± 

300 
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These measurements can be compared with results derived from 
Pioneer 10 Immersion and Pioneer 11 emersion at latitudes of 27° N and 
20° N, respectively. Woo and Yang [1978] concluded that p = 11/3 
between altitudes 500 and 3000 km, which is slightly larger than the 
values In Table 4.2, The discrepancy could be due to the different 
methods used for obtaining this parameter. The values in Table 4,2 are 
least-squares solutions for p from analysis of amplitude and phase 
spectra, whereas Wco and Yang [1978] presifctiably assigned a value on the 
basis of theoretical considerations and a visual inspection of measured 
amplitude spectra. In addition. Woo and Yang reached the same 
conclusions given above concerning the effect of collisions on the axial 
ratio. However, the VI results show that the axial ratio varies with 
altitude above 1200 km, which differs with the Pioneer observations of 
constant axial ratio in this region. Finally, Woo and Yang [1978] 
reported a fractional variation in electron density, 6N /N , of between 
9 and 18% near 1000 km altitude at 27°N. This result depends on the 
outer scale of the irregularities which they could not measure, but 
assumed to be 10 km. However, the Voyager phase measurements showed 
that the outer scale at 10° S was greater than 100 km, for which the 
Pioneer 10 measurement is scaled to a minimum of between 20 and 39% [Woo 
and Yang, 1976]; this is significantly larger than observed by VI. In 
this regard, the method adopted here for specifying 6N /N has the 
advantage of requiring no assumption about the outer scale of the 
irregularities, which is difficult to obtain. Furthermore, the VI 
observations extend the measurement of ^Ng/Ng to regions of anisotropic 
Irregularities (altitudes of 1200 to 3600 km), while the Pioneer results 
apply only for isotropic structures (near 1000 km altitude). 
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VI Emersion (0°, 314°M). The average electron density for this 
region of the Ionosphere was significantly lower than was measured by VI 
at Immersion [Eshleman et al., 1979b]. The spatial irregularities In 
electron density were too weak to produce appreciable scattering of the 
3.6-cm signal. However, weak scattering was observed at 13 cm, and 
amplitude spectra derived from the two smoothing intervals in Figure 4.5 
appear In Figure 4.6. The phase spectra for both smoothing intervals 
were dominated by noise. 



Figure 4.5 Profile of average electron concentration in the upper 
ionosphere at VI emersion. Regions 1 and 2 indicate altitude intervals 
for computing spectra [profile from Eshleman et al., 1979b]. 
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Figure 4.6 Amplitude spectra at 13 cm from the smoothing Intervals 
(regions 1-2) In Figure 4.5. See Figure 4.3 for additional comments. 


Table 4.3 Characteristics of Irregularities In Jupiter's upper 
Ionosphere at VI emersion. 

_3 

Height, Spectrum Y, Power-law Axial 6N , cm 

km type degrees Index ratio 

1500 13-cm a(t) 62 ± 4.6 3.3 ± 0.13 3.0 ± 0.6 1000 ± 200 

2500 13-cm a(t) 68 ± 6.1 3.0 ± 0.19 2.5 ± 0.6 590 ± 100 


Table 4.3 lists the least-squares solutions for parameters derived 
from these two amplitude spectra. The spectrum appears to folloi an 
inverse power law, with an exponent of about 3.0 to 3.3. Estimates of 
msjnetlc field direction, y, are mutually consistent. The fluctuations 
in electron density correspond to fractional variations ^Ng/Ng of 5 to 
10%. Again, the signal fluctuations in this frequency interval result 
from scattering by ionospheric Irregularities of size 2 to 100 
kilometers. 
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V2 Immersion (65**$, 260**M). The 13-cm signal experienced strong 
scattering as well as multiple-path propagation from steep gradients In 
electron density; no 13-cm spectra were reduced. The character of the 
fluctuations at 3.6 cm changed rapidly In this region; the two smoothing 
Intervals In Figure 4.7 Indicate the only segments during which the 
signal fluctuations appeared to be stationary. 



Figure 4.7 Profile of average electron concentration In the upper 
Ionosphere at V2 Immersion. Regions 1 and 2 indicate altitude intervals 
for computing spectra [profile from Eshleman et a1., 1979c]. 
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The amplitude and phase spectra at 3.6 cm from these Intervals 
appear In Figure 4.8, except that the phase spectrum of the first 
Interval corresponded to noise and Is not shown. Note that the 
amplitude spectrum obtained In the first Interval Is not characterlsitic 
of scattering from a power-law distribution of Irregularities. At high 
frequencies, the spectrum levels off at a value well above the noise, 
possibly as a result of scattering from the disturbances associated with 
an unstable plasma. In this regard. It Is perhaps significant that the 
Ionosphere at V2 Immersion lies near the magnetic footprint of the lo 
flux tube, where particle precipitation may affect the plasma dynamics 
[11. H. Acuna and N. F. Ness, private communication, 1979; Broadfoot et 
a1., 1979]. 

The residuals of the best theoretical fits to the spectra from the 
second smoothing interval were unexpectedly large. Apparently, the 
plasma Irregularities were not statistically homogeneous over the 
altitude range corresponding to the measurements, so the solutions for 
parameters derived from these spectra are not reported here. The 
limited amount of data available and the unique nature of Jupiter's 
ionosphere in this region restrict further Interpretation at this time. 
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Figure 4.8 Amplitude and phase spectra at 3,6 cm from the 
smoothing intervals (regions 1-2) in Figure 4.7. See Figure 4.3 for 
additional comments. No ionospheric parameters were obtained from these 
data. 
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V2 Emersion (50°S. 150°W). Two sets of 13-cm spectra, from the 
regions shown In Figure 4.9, appear In Figure 4.10. Figure 4.11 shows 
the squared coherency and cross spectnan phase at 13 cm from these 
regions. The signal fluctuations due to scattering at 3.6 cm exceeded 
measurement noise only for the amplitude spectrum of the lower altitude 
smoothing Interval. 



ELECTRON CONCENTRATION , cm'^ 


Figure 4.9 Profile of average electron concentration in the upper 
ionosphere at V2 emersion. Regions 1 and 2 indicate altitude intervals 
for computing spectra [profile from Eshleman et al., 1.979c]. 
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Figure 4.10 Amplitude and phase spectra at 3.6 and 13 cm from the 
smoothing intervals in Figure 4.9. See Figure 4,3 for additional 
comments. 
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Figure 4.11 Squared coherency and cross spectrun phase at 1? cm 
from regions 1 and 2 of Figure 4.9. See Figure 4.4 for additional 
comments. 


Table 4.4 Characteristics of irregularities In Jupiter's upper 
ionosphere at V2 emersion. 


-3 

Height, Spectrum Y, Power-law Axial 6N , cm 

km type degrees index ratio 


1500 

13-cm a(t) 

37 

t 

2.4 

3.5 

± 

0.05 

7.0 

± 

1.9 

3100 

± 600 

1500 

13-cm 4i(t) 

34 

± 

2.2 

3.5 

± 

0.04 

8.5 

± 

2.2 

4300 

± 

900 

1500 

3.6-cm a(t) 

30 

± 

1.9 

3.7 

t 

0.09 

8.0 

± 

2.1 

3600 

± 

700 

2500 

13-cm a(t) 

26 

± 

2.4 

3.4 

± 

0.07 

6.0 

± 

1.7 

1900 

* 

400 

2500 

13-cm 4>(t) 

28 

± 

3.2 

3.4 

± 

0.04 

8.0 

± 

2.3 

2500 

± 

bOO 
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Table 4.4 lists solutions for parameters derived from amplitude and 
phase spectra. The five estimates of y exhibit some spread but seem 
mutually consistent within the l-o uncertainties. The Inverse power-law 
exponent. p» ranges from 3.4 to 3.7. Spatial Irregularities of 
approximate size 2 to 140 kilometers contributed to scattering In the 
frequency range from 0.04 to 3 Hz. The spatial fluctuations In electron 
density correspond to fractional variations of 1 to 2% In Jupiter's high 
latitude Ionosphere. 

Solutions for ionospheric parameters Indicate a large value for the 
axial ratio over the full altitude range of appreciable signal 
scattering; this differs from the altitude dependence observed at the 
same heights in the equatorial ionosphere at VI Immersion (Table 4.2). 
The GSFC 04 model of Jupiter's magnetic field [Acuna and Ness, 1976] 
predicts a factor of nearly two between the field magnitudes In the 
ionosphere at V2 emersion (9.5 Gauss) and VI immersion (4.8 Gauss). 
Therefore, the magnetic field controls particle motion more effectively 
at high latitudes than at the equator, which may account for the strong 
anisotropy at both altitudes In Table 4.4, 

These measurements can be compared with results derived from 
Pioneer 10 emersion at latitude 59°N. Woo and Yang [1978] reported that 
p = 11/3 over the altitude range In Table 4.4. This Is again higher 
than values determined from the Voyager observations, although the 
comparison may be inappropriate since the measurements were made at 
similar latitudes but In opposite hemispheres (also see comments under 
VI immersion, above). The axial ratios in Table 4.4 agree with those 
measured during Pioneer 10 emersion. Woo and Yang [1976] deduced 
that ^Ng/Ng is 72% at 1500 km by assuming an outer scale of 100 km. 
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However, the V2 phase spectra Indicate an outer scale In excess of 140 
km, which would Increase the value of 6N„/N_ Inferred from the Pic. »r 
data. In either case, the Pioneer result Is considerably larger than 
values reported here, and seems unrealistic. Finally, Woo and Yang 
[1978] also suggested that the variation of Jupiter's magnetic field 
strength with latitude may lead to axial ratios that are larger at high 
latitudes than at the equator. The Voyager results contravene this 
conclusion -- although the field strength may affect the altitude 
dependence of the axial ratio, the largest axial ratio observed near the 
equator (10®S, Table 4.2} is similar to the largest axial ratio observed 
at higher latitudes (50“S, Table 4.4). 


CHAPTER 5 


SPATIAL IRREGULARITIES IN SATURN'S UPPER IONOSPHERE 
5,1 Introduction 

After their encounters with the Oovian system In 1979, Voyagers 1 
and 2 successfully obtained additional measurements of the Saturnian 
system In 1980 and 1981 [Stone and Miner, 1981, 1982]. This 
reconnaissance followed the first spacecraft flyby of Saturn by 
Pioneer 11 In 1979 [Dyer, 1980; Opp, 1980; Northrup et al., 1980]. 

Through analysis of dual -wavelength (13 and 3.6 cm) radio data from 
the occultatlons of Voyagers 1 and 2 by Saturn, Tyler et al. [1981a, 

1982] derived altitude profiles of the free electron density in Saturn's 
ionosphere, among other results. As was the case at Jupiter, rapid 
fluctuations were observed in the radio signals received from the 
occulted Voyager spacecraft; here, the scintillations are analyzed to 
determine morphology of plasma irregularities in Saturn's upper 
ionosphere. During the occultation of Pioneer 11 by Saturn, no effects 
due to scattering from the ionosphere could be detected within the 
limitations of the radio equipment [Woo and Armstrong, 1980]. 

The pf'xt section describes the experimental geometry and the 
Voyager occultation measurements; Section 5.3 gives the spatial 
characteristics of Saturn's ionospheric irregularities as deduced from 
the Voyager scintillation data. In Chapter 6, the magnetic field 
orientations at two locations are inferred from the alignment of the 
irregularities. These results have been published elsewhere [Hinson, 

1983] , 
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S.2 Experimental Geometry and Measurements 

Figure 5,1 shews the view from earth of the Voyager occultatlons by 
Saturn. Voyager 1 Immersion occurred on November 13, 1980, at 70°S, 
SOy^W, with solar zenith angle 6 = 91" (evening). The spacecraft 
emerged from behind the atmosphere at 2"S, 157"W, 6 =86" (early 
morning). The locations of immersion and emersion for the occultatlon 
of Voyager 2 on August 26, 1981, were 31°N, 36°W, 6 = 87" (evening) and 
26°S, 270"W, 6 o 93® (predawn), respectively. All longitudes above are 
expressed In the Saturn longitude system, or SLS [Desch and Kaiser, 
1981]. Radio-wave scattering from electron density structures 1n the 
upper Ionosphere caused scintillations In the signals received only 
during VI Immersion and V2 emersion; parameters which characterize these 
two occultatlons are listed In Table 5.1. Here, as before, v* 1s the 
ray path velocity in the plane of the sky; A indicates the orientation 
of V* (see Figure 6.1); and D Is the distance between the spacecraft and 
limb of Saturn. No scintillations due to scattering were observed 
during the occultatlons by the upper Ionosphere at VI emersion or V2 
immersion, probably as a result of the limited s1gna1-to-no1se ratios 
and the relatively low average electron densities at these locations; 
th«se data will not be discussed further. 
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VI 70°S 307'»W 17.2 142. Z.Zl 5.44 2.84 

V2 26“S 270"W 20.6 358. 2.07*10^ 5.21 2.72 

^P1 anetocentric 
^SLS CDesch and Kaiser, 1981] 



Figure 5,1 Schematic view from earth of the Voyager occultatlons 
by Saturn showing the paths followed by the spacecraft radio images. 



Saturn's Ionosphere 


It should be noted that scintillations were observed In connection 
with radio-wave propagation through the lower Ionosphere at all four 
locations. These relatively strong perturbations of the received signal 
are believed to be the result of multiple-path propagation through 
Ionized layers — similar effects were observed In the Voyager 
occultatlon data from Jupiter. In both cases, the weak-scattering 
theory of Appendix A and the models for Ionospheric irregularities In 
Chapter 2 are Inappropriate for Interpreting these data. The discussion 
In the remainder of this chapter pertains only to scattering from the 
topside Ionosphere. 

Again following the procedure for preliminary data reduction 
(Section 3.2), values of signal frequency, f(t), and amplitude, A(t), 
were measured each 0.064 s for VI Immersion at both wavelengths, 0.041 s 
for V2 emersion at 13 cm, and 0.027 s for V2 emersion at 3,6 cm. When 
displayed on an Argand diagram, the scintillations resulting from 
propagation through Saturn's Ionosphere are similar in appearance to 
those plotted in Figure 3.1 for VI immersion at Jupiter, with the 
exception that strong scattering was not observed at the 13-cm 
wavelength [cf, Hinson, 1933]. These values of f(t) and A(t) werp 
reduced by the methods In Sections 3.3 and 3.4 that were applied 
previously to data from Jupiter (Chapter 4). The temporal spectra for 
phase, <l)(t), and normalized amplitude, a(t), are discussed below, along 
with the solutions for the four ionospheric parameters. 
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5.3 Results 

VI Immersion (7Q”S, 307°W). Figure 5.3 displays two 13-cm 
amplitude spectra from the altitude intervals in Figure 5.2, with the 
best-fit theoretical spectra superimposed. The indicated noise levels 
were derived from measurements preceding ingress (see Section 3,4); 
weak scattering from ionospheric irregularities accounts for that 
portion of the scintillation spectra that exceeds the noise level. The 
13-cm phase spectra from both intervals, as well as all 3.6-cm spectra, 
were dominated by noise and ure not shown. 



Ficure 5.2 Average electron concentration as a function of 
altitude above the 1 bar pressure level for the upper ionosphere at VI 
immersion [profile from Tyler et al., 1981a]. Regions 1 and 2 Indicate 
altitude intervals for computing the spectra in Figure 5.3. 
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Least-squeres solutions for Ionospheric parameters were derived 
from the amplitude spectra; the results are listed In Table 5.il» with 
the 1-0 uncertainties. The altitude assigned to a measurement again 
corresponds to the center of the appropriate smoothing Interval. Above 
1 Hz. the decrease In the amplitude spectra Is consistent with an 
Inverse power-law distribution of Irregularities with exponent near 3.5, 
similar to values observed at Jupiter [see Section 4.3; Woo and Yang, 
1978] and the earth [Dyson et al., 1974; Phelps and Sagalyn, 1976]. As 
with Jupiter, Saturn's magnetic field orientation should not vary 
appreciably with radius over the extent of the upper Ionosphere, so the 
lack of variation In the estimates of y was expected. 

Recall that under the assumption of "frozen" Irregularities, there 
Is a direct proportionality between the frequency of signal fluctuations 
and the size of the scattering structures (see Section 4.3). Using this 
interpretation, Irregularities of approximate size 3 to 60 km contribute 
to scattering over the frequency range In Figure 5.3. These lower and 
upper limits result from noise limitations and the length of the 
measurement time interval, respectively, and do not represent either 
inner or outer scales. The spatial variations in electron density, 6N^, 
listed in Table 5.2 refer only to irregularities in this size interval, 
and correspond to fractional density variations of about 10%. 
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Amplitude 


noise leve 



V 

noise lev 

n 

Region 2 


FREQUENCY. Hz 


FREQUENCY, Hz 


Figure 5.3 Spectra of scintillations in 13-cm amplitude from the 
smoothing intervals in Figure 5,2. The l-o uncertainties apply at all 
sample frequencies. Theoretical spectra (smooth curves superimposed) 
were fitted in a least-squares sense to the measured spectra over the 
frequency range where the scintillations due to scattering exceeded the 
noise level by at least 3 dB. 


Table 5.2 Characteristics of irregularities in Saturn's upper 
ionosphere at VI immersion. 


Height, 

km 


Spectrum 
t 




decrees 


Power- law 
index 


13-cm 

a(t) 

29 + 2 

3.6 

+ 0.1 

13-cm 

a(t) 

27 ± 2 

3.3 

± 0.2 


Axial 

ratio 


6Ng, cm 


2600 + 500 
1300 ± 300 











Saturn's Ionosphere 


V 2 Emersion (ZS^S. 270*’W). Plasma Irregularities produ<;e<l weak 
scattering at both wavelengths during the occultation by the upper 
ionosphere. The amplitude and phase spectra appearing in Figure 5.5 
were computed over the altitude intervals in Figure 5.4; in the 
intervals where phase spectra are absent, noise dominated the 
measurements. The best-fit theoretical spectra are superimposed in 
Figure 5.5. 



ELECTRON CONCENTRATION, cm*^ 


Figure 5.4 Average electron concentration as a function of 
altitude above the 1 bar pressure level for the upper ionosphere at V2 
emersion [profile from Tyler et al., 1982]. Regions 1-3 indicate 
altitude intervals for computing the spectra in Figures 5,5 and 5.6. 
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For completeness and as a test of the theory, the squared coherency 
and cross spectrum phase from region 2 of Figure 5.4 are shown In Figure 
5.6. The high-frequency portion of the spectra corresponds to 
measurement noise for which the cross spectrum phase has a random 
distribution. Using the parameters obtained through analysis of the 
amplitude and phase spectra, the theoretical squared coherency and cross 
spectral phase were computed and superimposed In Figure 5.6. The 
parameter measuranents are consistent with the theory, as indicated by 
the general agreement between curves In each panel. 

Table 5.3 lists solutions for parameters derived from amplitude and 
phase spectra, with 1-a uncertainties. On the basis of the close 
agreement between the measurements and the theory, the spatial spectrum 
of Irregularities 1s consistent with an Inverse power law with exponent 
near 3.5. The four measurements of y are mutually consistent. These 
signal fluctuations were caused by anisotropic irregularities of 
approximate size 2 to 40 km (over the frequency Interval where the 
scintillations exceeded the noise). The values of 6N^ apply only for 
this range of sizes and correspond to fractional variations In density 
of about 20%. 

The measurements of axial ratio in Tables 5.2 and 5.3 show the same 
variation with altitude and latitude as was previously observed at 
Jupiter (see Tables 4.2 and 4.4). The axial ratio apparently increases 
with altitude at V2 emersion where Saturn's magnetic field Is relatively 
weak (0.24 Gauss), while the increased field strength at VI immersion 
(0.61 Gauss) Is believed to be responsible for the large axial ratios at 
both altitudes. (Both magnetic field magnitudes were calculated from a 
model proposed by Connerney et al. [1982]). 
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Figure 5.6 Squared coherency and cross spectrum phase at 13 cm 
from region 2 of Figure 5,4. The 1-a uncertainties are Indicated at one 
frequency. The high frequency portion of each spectrum corresponds to 
measurement noise for which the cross spectrum phase has a random 
distribution. Theoretical curves are superimposed at frequencies below 
4 Hz. 


Table 5.3 Characteristics of 1 rregularltles in Saturn's upper 
Ionosphere at V2 emersion. 


Height, 

km 

Spectrum 

type 

Y, 

degrees 

Power- law 
Index 

Axial 

ratio 

cm"^ 

2100. 

13-cm a(t) 

46 + 3 

3.6 + 0.2 

7 ± 3 

3500 + 700 

2100. 

13-cm <t>(t) 

43 t 4 

3.4 ± 0.2 

5 ± 2 

5500 ± 1100 

2100. 

3.6-cm a(t) 

41 ± 6 

3.6 ± 0.4 

4 + 2 

3600 ± 700 

2500. 

13-cm a(t) 

42 ± 3 

3.5 ± 0,4 

11 ± 6 

1200 ± 300 
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CHAPTER 6 


MAGNETIC FIELD ORIENTATIONS IN THE IONOSPHERES OF JUPITER AND SATURN 

6.1 Introduction 

In this chapter, magnetic field orientations In the Ionospheres of 
Jupiter and Saturn are Inferred from the occultatlon measurements In 
Chapters 4 and 5. Section 6.2 describes the procedure for deducing the 
field direction; Sections 6.3 and 6.4 report results for Jupiter and 
Saturn, respectively. Comparisons with the predictions of magnetic 
field models are Included for both planets. The chapter ends with a 
discussion of results In Section 6.5. 

In summary, at Jupiter, the occultatlon measurements conflict with 
the predictions of current field models, which apparently do not 
accurately represent the field near the planet. However, when used as a 
test of Saturnian magnetic field models, the occultatlon measurements 
generally confirm the model predictions, but also suggest that a small 
adjustment of model parameters might Improve the field representation. 

6.2 The Inference of Magnetic Field Orientation from Occultatlon 
Measurements 

It has been assumed here that the Ionospheric Irregularities are 
field aligned and symmetrical with respect to rotations about the field 
direction (cf. Section 2.3.1). Hence, the measurements of the 
orientation, y, 1n Tables 4,2-4. 4, 5.2, and 5.3 are related to the 
magnetic field directions In the Ionospheres of Jupiter and Saturn. To 
understand the constraint provided by these measurements, recall that 
the scintillation spectra are sensitive only to refractive variations in 
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the plane that Is perpendicular to the direction of microwave 
propagation; the spectra are essentially unaffected by variations along 
t,he propagation path (cf. Appendix A). Consequently, only the field 
component that Is perpendicular to the direction of propagation will 
Introduce anisotropy into the plasma structures In a way that can 
Influence these occultatlon measurements. The component of the 
planetary field (In a spherical r, 6, 0 representation) was nearly 
tangent to the propagation path for each of the Voyager occultatlons by 
Jupiter and Saturn. This geometry allows a determination of the 
relative sizes of the r and 6 field components, but no information 
concerning the <(i component can be obtained through scintillation 
techniques. Although the componeni affects the angle cf polarization 
of the microwave signals through the phenomenon of Faraday rotation, the 
right-hand circular polarization of the signals transmitted by Voyager 
was unsuitable for observations of this effect. 

There are two ambiguities In the reduction of these observations to 
field orientation. First, a reversal of field direction would not 
affect the ionospheric irregularities; this results in a 180° 
uncertainty in y. Second, the theoretical spectra are sensitive only to 
the absolute value of y. which allows two possible mirror image 
solutions for the orientation [cf. equations A36 and 2.4]. Both 
ambiguities were resolved by choosing the field orientation that 
minimizes the difference between the occultatlon measurements and 
current Jovian or Saturnian magnetic field models. Other possible 
solutions for the orientation are well separated from the model 
predictions, which can be seen by considering Figure 6.1, below, and the 
angles in Tables 4.1-4. 4, 5, 1-5.3, and 6. 1-6. 2. 
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To facilitate comparisons with magnetic field models In the next 
two sections, the field orientations Inferred from occultatlon 
measurements are expressed In terms of r, the angle measured clockwise 
relative to the outward radial direction on the visible limb of the 
planet, as defined in Figure 6.1. This angle follows directly from the 
measurements of y (Tables 4. 2-4. 4, 5.2, and 5.3) through the angle 4, 
which defines the geometry of the occultations (Tables 4.1 and 5.1). At 
Saturn, a small additional correction to r (-1® for VI, +2® for VZ) has 
been Included to account for the motion of the Ionosphere across the 
propagation path that results from Saturn's rotation and the slight 
Inclination of the rotation axis at the time of the occultations. Due 
to the small inclination of Jupiter's rotation axis during the Voyager 
encounters, this correction Is negligible at Jupiter. 
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FigtJre 6.1 View from earth of the visible limb of the occulting 
planet (Jupiter or Saturn) showing the raypath velocity, v*, and the 
component of the planetary magnetic field, that is perpendicular to 
the line of sight. The angles A and r are measured clockwise relative 
to the outward radial direetior, r, at the occultation point. 
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6.3 Magnetic Field Orientations In Jupiter's Ionosphere 

Although numerous models for Jupiter's magnetic field have been 
proposed [cf. Smith et al., 1976], only four will be considered here; 
each Is based on magnetometer measurements obtained during the Pioneer 
11 (Pll) flyby. Acuna and Ness [1976] fitted a third-order spherical 
harmonic expansion to the flux gate magnetometer (FGM) data recorded 
between 1.7 and 6 Jupiter radii to produce Model 1 (GSFC 04). In this 
representation, the observed magnetic field Is attributed entirely to 
sources Internal to Jupiter; external sources are neglected. Smith et 
al. [1976] derived Model 2 (Pll 31 IE) through analysis of the Pll 
vector helium magnetometer (VHM) data. This spherical harmonic 
expansion includes a dipole repr seitatlon for the sources external to 
Jupiter In addition to a third-order expansion for Internal sources. 
Connerney [1981] proposed Models 3 and 4 after Independently analyzing 
the Pll FGM and VHM data, respectively. In the following manner. First, 
the Pioneer 10, VI, and V2 magnetometer data were used to construct a 
model for the magnetic field arising from currents In the Jovian 
magnetosphere. Next, this model was applied to remove the contribution 
by external sources from the Pll data. When used as Input for a 
generalized Inverse analysis, the resultant field yielded a third-order 
spherical harmonic model of Jupiter's Internal field for each data set. 

Table 6.1 summarizes the results for the magnetic field orientation 
in Jupiter's Ionosphere, and Includes comparisons with the four models 
just described. The third column lists occultation measurements of r at 
the locations specified In the first two columns. Columns 4 to 7 give 
predictions of r from the four Jovian magnetic field models. Relative 
to the 1-0 measurement uncertainties, the differences between the 
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occultatlon measurements and the predictions of all four models are 
significant. No discussion of the Pioneer occultatlon measurements of 
field orientation [Woo and Yang, 1978] Is Included here, since those 
results did not Include measurement uncertainties, and comparisons with 
magnetic field models are therefore Inconclusive. 

As mentioned previously (Section 4.3), the field orientation should 
not vary appreciably with radius over the vertical extent of Jupiter's 
upper Ionosphere. Nevertheless, in response to the longitudinal 
structure In the magnetic field, the orientation does vary slightly with 
position along the radio path within the ionosphere. This path 
dependence has been Included In computing an effective r from the four 
models; these values and the measurements in Table 6,1 should be 
Interpreted as local spatial averages of the field orientation In 
Jupiter's ionosphere. However, the total variation along the raypath 
is not larger than a few degrees, this effect cannot account for the 
differences between the predictions of the field models and the 
occultatlon measurements. 
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Table 6.1 Jupiter's magnetic field orientation at ionospheric heights. 

r , Magnetic Field Angle, degrees 

OccuVtatlor Connerney [1981;] 


Lat.^ 

Lonq.^ 

Measurement 

GSFC 04 

Pll 31 IE 

FGM 

VHM 

10"S 

63°W 

119 

i 

2.4 

136 

136 

137 

132 



119 

± 

1.8 







117 

± 

1.7 







117 

± 

1.7 







119 

± 

1.0 







117 

± 

1.6 





0® 

314°W 

303 

± 

4.6 

284 

288 

283 

290 



297 

± 

6.1 





50“S 

150®W 

200 

i 

2.4 

196 

196 

197 

197 


203 + 2.2 
207 ± 1.9 
211 ± ?.4 
209 ± 3.2 

^Planetocentric 

^Sys III, 1965 rseidelmann and Divine, 1977] 
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6.4 Magnetic Field Orientations In Saturn* s Ionosphere 

The first models of Saturn's magnetic field were derived from Pll 
magnetometer data, and consisted of a dipole aligned with the rotation 
axis but with Its center offset northward by 0.04 R^, where R^ Is 
Saturn's radius [Smith et al., 1980a, b; Acuna and Ness, 1980; Acuna et 
al., 1980]. Spacecraft observations of the charged particle absorption 
signatures of Saturn's moons and rings also supported the "offset 
dipole" (OD) model [Chenette and Davis, 1982]. Through subsequent, 
analysis of the combined VI and V2 magnetometer data, Coniierney et al. 
[1982] derived the "Z3" model for Saturn's field, which consists of the 
first three zonally symmetrical terms In a spherical harmonic 
expansion. Aside from a nonzero octupole coefficient, g^®, in the Z3 
expansion, the two models are essentially the same -- both are 
symmetrical about Saturn's rotation axis and include no 4» component. As 
the Importance of the octupole term decreases rapidly with Increasing 
distance from Saturn, the charged particle absorption signatures are 
also consistent with the Z3 model [Acuna et al., 1982]. 

Table 6,2 summarizes results for the magnetic field orientation in 
Saturn's ionosphere at the locations specified in the first two 
columns. The third column lists occultation measurements of r, the 
angle of Saturn's field component perpiL'dicular to the line of sight as 
defined in Figure 6,1, For comparison, predictions of r from the OD and 
Z3 magnetic field models appear in the last two columns. The difference 
between the two theoretical values at each location illustrates the 
effect of gj° on the field representation at ionospheric heights. At 
both locations in Table 6.2, the occultation measurements are 
intermediate to the predictions of the two models; within the 
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mecsurement uncertainties, each model is in reasonable agreement with 
the occultation results (also see Figures E3 and E4 of Appendix E). 

As mentioned previously, the occultation measurements in Table 6.2 
should be regarded as local spatial averages. The systematic variation 
of r along the radio path in Saturn's ionosphere was calculated on the 
basis of both the 00 and 23 models and the geometry of the Voyager 
occultations. Unlike the case at Jupiter, this effect was found to be 
insignificant when compared with the measurement uncertainties. 


~^able 6.2 Saturn's magnetic field orientation at ionospheric heights. 

r, Magnetic Field Angle, degrees 
, „ Occul tat i on 


Lat.^ 

lor.g.^ 

Measurement 

23 

00 

70“S 

307“W 

170 t 2 

166 

171 



168 ± 2 



26 "S 

270“W 

226 t 3 

229 

221 


223 ± 4 

221 t 6 

222 ± 3 


^Planetocentric 

^SLS [Desch and Kaiser, 1981] 
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6.5 Discussion 

When comparing the occultation observations to the model 
predictions at Jupiter (Table 6.1) and Saturn (Table 6.2), It Is 
Important to note certain limitations of the spacecraft magnetometer 
data [Connerney, 1981]. First, at Jupiter, the Pll FGM and VHM data are 
incompatible in the sense that systematic differences between the two 
sets of field measurements are unexpectedly large. Second, when 
deriving field models from magnetometer data, the representation for the 
internal field is sensitive to assumptions about the fields of external 
origin and to the alternative methods of data inversion; consequently, 
there are substantial differences among the four Jovian field models 
mentioned in Section 6.2 (see Figures 14 and 15 of Connerney [1981]). 
Finally, magnetometer measurements along a spacecraft flyby trajectory 
are insufficient to produce a unique model of a planet's internal 
magnetic field. Because of these limitations, complementary 
measurements, such as the field orientations inferred from Voyager radio 
occultations, could provide valuable information on the Jovian and 
Saturnian fields. 

Jupiter's magnetic field is more complex than that of a simple 

dipole — it has significant quadrupole and octupole moments [Smith and 

Gulkis, 1979]; the Same may be true of Saturn's magnetic field 

[Connerney et al., 1982]. In a spherical harmonic representation, 

-3 -4 

dipole terms vary as r , quadrupole terms vary as r , and octupole 

-5 

terms are proportional to r . Consequently, the quadrupole and 
octupole terms in the expansion are relatively more important in the 
ionosphere than at the remote position of spacecraft magnetometer 
measurements. Therefore, in principle, the field orientations inferred 
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from radio occultatlons can provide a sensitive test for deviations from 
zonal symmetry In the magnetic field, and a basis for Improving the 
higher-order coefficients in field models. 

Unfortunately, occultatlon measurements at three locations for 
Jupiter and two locations for Saturn provide only a weak constraint on 
the twelve quadrupole and octupole terms In the respective fields. The 
difficulties In modifying Saturnian field models are compounded by the 
current uncertainty In Saturn's rotation rate, which can Introduce t;p to 
±40° of error In longitude during the Interval between the VI and V2 
encounters [Desch and Kaiser, 1981], Furthermore, in adjusting the 
model parameters, the residuals of both the occultatlon measurements and 
magnetometer data must be considered, which is beyond the scope of this 
research; no specific revisions to present field models are proposed 
here. 

Nevertheless, two general conclusions can be drawn from the results 
in Tables 6.1 and 6.2. First, at Jupiter, current models are apparently 
Inadequate f'?r accurately representing the field near the planet. The 
large discrepancies between the occultatlon measurements and model 
predictions in Table 6.1 may be the result of complicated but as yet 
unknown fine structure In Jupiter's magnetic field. Second, at Saturn, 
a small adjustment of the coefficients In current field representations 
would remove the differences between the model predictions and 
occultatlon measurements. For example, a field model with a 

coefficient Intermediate to the values in the Z3 (2743 nanotesla) 
and OD (0 nanotesla) representations would agree more closely with the 
occultatlon results at both locations in Table 6.2. A 20% reduction of 
in the Z3 model would change the predictions for r by +1° at 70°S 
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and -2® at 26*S. Alternatively, the 93 ° coefficient In the Z3 model 
could be retained, a possibility which Is preferred by Northrop and Hli?, 
[1983], and the discrepancies In Table 6.2 could be reduced by assigning 
small nonzero values to other octupole or higher order coefficients 
94 °) in the field expansion. 
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INTERNAL GRAVITY WAVES IN TITAN’S ATMOSPHERE 
7.1 Introduction 

The flyby of the Saturnian system by Voyager 1 Included a radio 
occultation by Titan on November 12, 1980 [Stone and Miner, 1981]. 

Tyler et al. [1981a] and Lindal et al. [1983] analyzed the average 
effect of Titan's atmosphere on the two radio signals transmitted by the 
spacecraft and received on the earth, and thereby derived vertical 
profiles of gas refractivlty , molecular number density, pressure, and 
temperature. Through subsequent analysis of the occultation signals, 
Lindal et al. [1983] also showed that there was neither measurable 
microwave absorption In the neutral atmosphere nor a detectable free 
electron density In the Ionosphere. This chapter considers the 
relationship between rapid signal fluctuations during the occultation 
and small-scale Irregularities In Titan’s neutral atmosphere. These 
scintillation data provide the only Information presently available 
concerning atmospheric structure and dynamics on length scales from a 
few hundred meters to a few kilometers In the troposphere and 
stratosphere. 

Section 7.2 describes the experimental geometry of the Voyager 
occultation by Titan and the fundmental measurements. In Section 7.3, 
spatial characteristics of the atmospheric Irregularities are determined 
from the reduced scintillation data; In addition, the altitude 
dependence of the intensity variance Is calculated. Section 7.4 
interprets these results, and confirms the hypothesis that freely 
propagating Internal gravity waves are present In Titan's upper 
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atmosphere. In addition, wave parameters, Including the vertical fluxes 
of energy and momentum, are estimated, and the source of the gravity 
waves Is considered. These results have been published elsewhere 
[Hinsor and Tyler, 1983]. 

7.2 Experimental Geometry and Measurements 

Figure 7,1 shows a trajectory pole view of the Voyager 1 
occultatlon by Titan; the direction to earth Is upward. Occultation 
immersion occurred near latitude 6®N, solar zenith angle 6 = 94® 
(evening), with emersion about 11 minutes later at 9“S, 6 ° 86® (early 
morning). (The orientation of Titan's pole Is currently unknown, but is 
assumed to be normal to the mean orbital plane [Davies et al., 1980]. 
Hence, these latitudes are approximate,) Before 1mjr.ers1on, the 
spacecraft antenna boresight was offset 0,11° In order to compensate for 
this amount of refraction. This configuration allowed complete 
amplitude and frequency measurements at both radio wavelengths for the 
tenuous upper part of Titan's atmosphere. Before emersion, the antenna 
boresight was repositioned to be 2.36° from ear.h line, which then 
allowed amplitude and frequency measurements at the 13>cm wavelength 
over the maximum possible range In pressure. However, the 3.6-cm data 
were of limited value In this orientation because of the relatively 
narrow antenna beamwidth at this wavelength. (Half-power width of the 
spacecraft antenna pattern is 2.5° at 13 cm and 0.6° at 3.6 cm). For 
this reason, the 3.6-cm data from emersion will not be discussed 
further; however, these data should be considered In future work, as the 
signal was detected briefly. These antenna positions were chosen to 
Insure good results despite significant a priori uncertainties 1n the 
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size, orbital position, and atmospheric density of Titan [Tyler et al., 
1981a; Lindal et al., 1983]. 

In order to be consistent with the notation used In the literature, 
the observations during the occultatlon by Titan's neutral atmosphere 
are reported In terms of signal intensity rather than amplitude (cf. 
Section 3.3). As discussed In Section A. 7 of Appendix A, this change Is 
Insignificant for the case of weak scattering. 



Figure 7.1 Trajectory pole view of Titan. Direction to earth Is 
upward and the antenna orientations for immersion and emersion are 
indicated. 
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Values of signal frequency, f(t), and Intensity, I(t), were 
measured (Section 3.2) each 0.041 s for Immersion at 13 cm, 0.027 s for 
Immersion at 3.6 cm, and 0.082 s for emersion at 13 cm. Limitations of 
the radio equipment and observation time resulted In standard deviations 
In the estimates of f(t) and I(t) of approximately 0.7 Hz and 7%, 
respectively, before Immersion, and 1.5 Hz and 23% after emersion. 

These values of I(t) are used below to plot the time history of 
Intensity (Figure 7.2) and to calculate scintillation spectra (Figure 
7.3). The signal parameters also were measured at one half this 
sampling rate. These alternate values of I(t) have reduced noise 
effects, and therefore are preferable for determining the altitude 
profile of the Intensity variance (Figure 7.5). 

Intensity profiles for Immersion at 3.6 and 13 cm, and emersion at 
13 cm appear In Figure 7.2. To facilitate comparison, the time scale 
for emersion has been reversed and compressed to compensate for 
geometrical differences between Immersion and emersion. In the first 
segment of each profile, where the average Intensity remained constant, 
the radio path from spacecraft to earth was above Titan's atmosphere. 
During the occultation, the general shape of each curve was determined 
by the combined effects of atmospheric refraction and antenna pointing 
[cf. Lindal et a1., 1983; also see Appendix D]. The 13-cm signal was 
received throughout both immersion and emersion except when It was 
blocked by Titan's surface, while the 3.6-cm signal was lost In the 
noise because of the antenna pattern effects before the surface was 
reached. The altitude of ray perlapsis relative to Titan's surface Is 
indicated in Figure 7.2 and applies to all three profiles [private 
communication, G. F. Lindal, 1982]. The last portion of each curve 
corresponds to the measurement noise level. 
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Figure 7.2 Profiles of Intensity vs. time showing Immersion for 
3.6 and 13 cm, and emersion for 13 cm. The time scale for emersion has 
been reversed and compressed to account for geometrical differences 
between immersion and emersion. The combined effects of atmospheric 
refraction and antenna pointing determined the average shape of each 
curve. The 13-cm signal was received throughout immersion and emersion 
except when It was obscured by the surface; the 3.6-cm signal was lost 
in the noise because of the antenna orientation before the surface was 
reached. Titan's neutral atmosphere caused Intensity fluctuations on 
time scales from about 0.1 to 1.0 s at both wavelengths, indicating the 
presence of variations In refractivity on length scales from a few 
hundred meters to a few kilometers at all altitudes within about SO km 
of the surface. 
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During emersion, the antenna offset closely coincided with the 
angle of refraction by Titan's lower atmosphere. In this case, the 
signal extinction by the surface Is clearly evident. The cut-off by the 
surface during Immersion Is unambiguous at a lower sampling rate which 
Increases the signal -to-nolse ratio [Lindal et a1., 1983]. 

Small fluctuations In intensity appear In the first segment of each 
of the curves In Figure 7.2. These rapid variations result from thermal 
noise, and are representative of the measurement uncertainties In the 
upper atmosphere. Titan's neutral atmosphere caused Intensity 
fluctuations which exceed the noise on time scales from about 0.1 to 
1.0 s at both wavelengths. As discussed In detail below, this Indicates 
the presence of variations In refractivlty on length scales from a few 
hundred meters to a few kilometers at all altitudes within about 90 km 
of the surface. The strongest scintillations were observed with the 
lowest point on the radio path between about 5 and 25 km above the 
surface. Weaker intensity scintillations occurred with ray perlapsis 
between the surface and 5 km, and between 25 and 90 km. Note that the 
scintillations show similar variation with altitude for Immersion and 
emersion — evidence that the atmospheric structures that cause these 
effects are not a localized phenomenon, as the observations were on 
opposite sides of Titan. 

The Intensity data were reduced through a two step procedure 
(Section 3,3). First, the technique of spectral analysis was applied to 
the weak Intensity fluctuations caused by scattering from the upper 
atmosphere. Second, the fractional variance of Intensity, (see 
equation 3.4), was calculated for data obtained throughout the 
occultatlon. These results appear in the next section. Analysis of 
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phase fluctuations was restricted by the slgnal-to-noise ratio of the 
data and the limitations of the weak-scattering theory. 

In order to Interpret the Intensity spectra, theoretical spectra, 
w^, were calculated from equation A41 In Appendix A (as discussed In 
Sect1''jn A. 7, w- ° 4w for weak scattering). It Is Important to note 

* A 

that the parameters 1> and <t> (see Appendix D) In this equation affect 

y ^ 

the shape of w^ [Haugstad, 19793, and vary with time during the 
occultation. Their values were determined using the formulas In 
Appendix D and parameters supplied by G. F. Lindal [private 
communication, 1982]. For all data considered here, 1 so that limb 
curvature had a negligible effect. However, <J)^ varied from 1.0 before 
immersion to less than 0.1 with ray periapsis near the surface, and the 
steady variation of during the time Interval of one spectrum had the 
effect of smoothing the measurements over a fractional range 1n 
frequency of about ±10%. To account for this frequency smearing, four 
theoretical spectra were calculated from equation A41 using four 
slightly different values of Tnen, these spectra were averaged to 
produce a smoothed theory, w^ , for comparison with the measurements. 
Tyolcal values of parameters used in these computations were 
V = 0.1 km s"^, v = 7.98 km s“^, = 1.01, 4» = 0.29, and 

J ^ J ^ 

D = 9150 km (see Appendix D for definitions); these values apply for 
scattering from Titan's atmosphere near 44 km during immersion. 

Comparisons of measured spectra with the smoothed theory, w^ . 
produced solutions for four parameters that characterize the scattering 
Irregularities In the neutral atmosphere (see Section 3,4). In this 
case, the parameters are (Section 2.3.2): (1) the axial ratio, e/f; 

(2) the orientation, (3) the Inverse power-law exponent, p, of the 
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spatial spectrum; and (4) the magnitude of the spatial variations In 
refractive Index, n^. The covariance matrix of the solution for the 
four parameters was computed by propagation of errors {Section 3.4) 
using the theoretical expression for as a model for the process. The 
parameter measurements and the uncertainties are dicussed In the next 
section. 

7.3 Results 

intensity spectra at 3.6 and 13 cm appear In Figure 7.3. Figure 

7.4 shows the altitude Interval of these measurements, which correspond 

2 2 

to weak scattering (a^ = 0.25 at 3.6 cm and = 0.076 at 13 cm; cf. 

Appendix B). The noise levels were determined from observations before 
Immersion (see Section 3.4) with a correction for the decrease In 
average signal Intensity. As this noise level was subtracted before 
fitting, large variations appear in the high frequency part of each 
measured spectrum. The best fit theoretical spectra are given as the 
smooth curves In Figure 7.3. When compared to the noise level, the 
fluctuations In 13-cm Intensity from- above 25 km during emersion were 
too small for useful spectral analysis. 
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figure 7.3 Weak-scattering intensity spectra ^t 13 cm (upper 
pane'i) and 3.6 cm (lower panel) during immersion from the smoothing 
interval centered at 44 km in Figure 7.4. Theoretical spectra (smooth 
curves superimposed) were fitted in a least-squares sense to measured 
spectra over the frequency range where the intensity fluctuations exceed 
the noise level by at least 3 dB. Note that the noise level was 
subtracted before fitting, resulting in large variations in the high 
frequency portion of both spectra. 
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TEMPERATURE, K 

Figure 7.4 Temperature-pressure profile with altitude scale for 
Titan's atmosphere at immersion [Lindal et al., 1983], Vertical bar 
indicates the smoothing interval for computiri the intensity spectra in 
Figure 7.3. 


Tab^ p 7.1 Characteristics of irreglarities in Titan's upper atmosohere 
at liinnersion. 


Height, km 
44 
44 


cm degrees 


Axial 

ratio 




-4 


3.6 

13. 


0. i 10. 
10. ± 20. 


3.7 ± 0.3 
3.1 ± 0.4 


4 ± 2 
6 ± 4 


(8 ± 3) • 10 
(9 ± 4) • 10 
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The solutions for atmospheric parameter? derived from the Intensity 
spectra are listed In Table 7.1. Each measurement should be Interpreted 
as a spatial average over the smoothing Interval In Figure 7.4. The 
uncertainties, which correspond to l-a, result from noise limitations 
and the steady variation of over the time period of the measurements, 
which tends to suppress the distinguishing features In the spectra. 

Recall that e/f and 4» characterize the shape and orientation, 
respectively, of the atmospheric Irregularities (see Section 2.3.2 and 
Figure 2.3). Applying this Interpretation to the present case 
(e/f ~4, (J< ~0), both horizontal and vertical variations In 
refractivlty probably were present In Titan's upper atmosphere. The 
structures were not Isotropic, but Instead had a tendency toward 
horizontal stratification. In addition, the normalized variations In 

_3 

refractive Index, n,./n„, were of order 10 . 

t o 

The transverse raypath velocity at 44 km altitude Is given by 

jv^i « 2.3 km s”^ (see Appendix D and the parameter values given 

above). This value greatly exceeds both the phace velocity of 

atmospheric waves and the bulk velocity of fluid motion which are likely 

to be encountered in Titan's atmosphere (see below). Hence, the 

Intensity fluctuations In Figure 7.3 result primarily from motion of the 

radio path across a fixed field of refractive Irregularities. Using 

this Interpretation, the scale sizes In the atmosphere which contributed 

to the spectral range 1 to 10 Hz lie between ajout 0.2 and 2 km. 

2 

The variance of normalized Intensity, , was measured over time 
Intervals of 1.3 s for Immersion at 13 cm, 0.87 s for Immersion at 

3.6 cm, and 2.6 s for emersion at 13 cm; the three curves appear In 

2 

Figure 7.5. In general, larger values of were observed at lower 
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2 

altitudes. Hence, each value of was assigned the minimum periapsis 
altitude in the corresponding measurement Interval. As mentioned 
previously, and v* varied smoothly as a function of periapsis 
altitude during the occultation. Each parameter was reduced to about 
10% of its free-space value when the raypath reached Titan's surface. 
Consequently, the resolution of the measurements improved from about 
8 km at 100 km altitude to about 0.8 km at the surface during immersion, 
with roughly half the resolution at the same altitudes during 
emersion. The standard deviation of was calculated by established 
methods [Jenkins and Watts, 1968], and appears at two altitudes in each 

curve. The smooth curves show the noise levels, which vary as 

-1 2 
<I(t)> ". At altitudes below about 90 km, clearly exceeded the 

noise — the observed scintillations resulted primarily from scattering 

fin Titan's atmosphere (cf. Figure 7.2). 

The measurements in Figure 7.5 quantitatively illustrate the 

general features mentioned previously in describing Figure 7.2. During 

immersion, spatial irregularities in Titan's atmosphere caused week 

Intensity scintillations when ray periapsis was between the surface and 

2 

about 5 km, or between 25 and 75 km, where was nearly constant. 
Strong scintillations occurred with ray periapsis between 5 and 25 km. 
(The basis for distinguishing between weak and strong scintillation is 
given in Appendix B.) At both wavelengths, increased by more than 
an order of magnitude when ray periapsis descended from 90 to 70 km. 
During emersion, o/ showed the same general characteristics as for 
immersion, although the increased noise level restricts 
interpretation. The theoretical curves, shown superimposed, are 
explained in the next section. 
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VOYAGER 1 /TITAN 



Figure 7.5a Normalized intensity variance at 13 cm (heavy line) as 
a function of perlapsis altitude during Voyager 1 Immersion at T' tan. 

The noise level varies as <I(t)>"^ and is indicated by the smooth 
curve. Altitude resolution improved from about 8 km at 100 km altitude 
to about 0.8 km at the surface. Uncertainties are shown at two 
altitudes. The theory (superimposed light line) is based on a model for 
the propagation of internal gravity waves through the vertical 
stratification in Titan's atmosphere. 
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VOYAGER 1 /TITAN 



Figure 7.5b Normalized intensity variance at 3.6 cm during 
Voyager 1 immersion at Titan. See Figure 7.5a for additional comments. 


VOYAGER 1 /TITAN 



Figure 7.5c Normalized intensity variance at 13 cm during 
Voyager 1 emersion at Titan. See Figure 7.5a for additional comments. 
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7 .4 Discussion 

Upper Atmosphere (above 25 km). The parameters In Table 7.1 

provide evidence for horizontal variations In refracti vlty , and hence 

dynamics, near 44 km In Titan's atmosphere. Above 25 km, the altitude 
2 

profile of (Figure 7,5) closely agrees with the predictions of a 
simple theory based on the chracterlstlcs of freely propagating 
atmospheric gravity waves [Hines, I960; Gossard and Hooke, 1975], as 
will be shown below. On the basis of these observations, the 
scintillations that occurred during occultatlon by the upper atmosphere 
are primarily the result of scattering from the density disturbances 
associated with gravity waves. (Note that the high vertical stability 
of the atmosphere above 25 km excludes convective turbulence as an 
explanation for Intensity scintillations.) 

To understand the altitude variation of , consider a parcel of 
gas displaced adlabatical 1y In the vertical direction. During a gradual 
displacement, the fractional difference In density between the displaced 
parcel and Its new surroundings 1s given by [cf. Gossard and Hooke, 

1975, p. 71-75]: 


op _ 1 56 

p " 0 5z 


dz 


(7.1) 


where dp Is the density difference, p 1s the density of the undisturbed 
atmosphere, 0 1s the potential temperature, and the vertical 
displacement, dz. Is assumed to be small compared to the atmosphere's 
scale height. 
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The derivative of 6 can be expressed as [cf. Turner, 1973, p,5]: 


1 d6 1 , dT 

“S" ~5T 


-3-J 5 n 


(7.2) 


where T Is the temperature of the undisturbed atmosphere, g Is the 
gravitational acceleration, and c^ Is the specific heat per unit mass at 
constant pressure. The density variations, 6p, result from the combined 
effects of adiabatic expansion and the vertical stratification of the 
undisturbed atmosphere. 

To approximate the behavior of atmospheric gravity waves, It Is 

assumed here that 6z Increases exponentially with altitude In proportion 

to '50 as to maintain a constant vertical energy flux [Turner, 

1973, p. 10-11, 22-23]. As discussed In Appendix A (see Section A. 7 and 

2 

the paragraph following equation A40), Is proportional to both the 

2 

variance of refractive Index, <n^ >, and the variance of fluctuations in 
atmos; heric density. By combining equations 7.1 and 7.2 the theoretical 
intensity variance in this model varies with altitude as 



(7.3) 


The constant of proportionality depends on experimental geometry, radio 
wavelength, and the vertical energy flux and spatial characteristics of 
the wave. 

2 

The theoretical was calculated from equation 7.3 using values 
of p, T, and &T/dz measured during the radio occultatlon [Lindal et al., 
1983], values of g{z) determined from Titan's mass [Tyler et al., 

1981a], and a value of Cp appropriate for Ng. The theoretical curves 
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are superimposed In Figure 7.5, The constant of proportionality In 
equation 7.3, which amounts tc an additive offset on the logarithmic 
plot, was adjusted to obtain a good fit. The ratio of these constants 
at 3.6 and 13 cm during Immersion, about 3, Is consistent with 
scattering from a power-law distribution of Irregularities with exponent 
p • 4.3 (see Section A. 7 and the paragraph following equation A40). 

This value Is within 3o of the estimates In Table 7.1. 

Above 25 km during immersion, the altitude variation of the 

observed scintillations Is entirely consistent with a hypothesis of 

freely propagating gravity waves; Interpretation for emersion is 

somewhat restricted by noise, but also appears to be consistent with 

this hypothesis. Note the detailed agreement between the theory and 

3.6-cm measurements (Figure 7.5b): both exhibit a nearly constant level 

of scintillation between 25 and 60 km, a local maximum at 65 km, and a 

rapid decrease above 70 km that Includes several changes In slope. A 

region of high vertical stability, where dT/az Is large (Figure 7.4), 

2 

leads to the enhancement of near 65 km, while the abrupt decrease 
above 70 km results from the atmosphere's average vertical structure, 
not from wave attenuation or reflection. 

The model described by equations 7.1 to 7.3, though not rigorous, 
has the advantage of requiring relatively few assumptions while still 
Incorporating the principal features of gravity-wave propagation. In 
this regard, the level of analysis attempted here Is consistent with the 
limitations of the available data. To avoid unnecessary complications 
In the discussion below. It Is assumed that a single coherent gravity 
wave caused the scintillations. In reality, the signal fluctuations are 
probably the result of scattering from a superposition of gravity waves. 
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since the spectra In Figure 7.3 reveal a continuous distribution of 
energy over a range of length scales. Nevertheless, this assumption 
should be acceptable for order-of*magn1tude calculatlortiy. 

The physical characteristics of the gravity waves can be Inferred 

from the parameters listed In Table 7.1. Near 44 km, n = 2»10"^ m"^, 

-3 

p => 0.5 Kg m [Lindal et al,, 1983], and the fractional variations In 

-3 

density associated with the wave must satisfy 6p/p = n^/n^ « 10 (for 

the relationship between refractlvlty and density, see Section 2.2). 

The horizontal and vertical wavelengths are roughly 4 and 1 km, 

respectively, which correspond to the observed size and axial ratio of 

the Irregularities. Then, using standard formulas for the linearized 

relationships among field variables for gravity waves [Gossard and 

Hooke, 1975, p. 75], the horizontal and vertical fluid velocities 

,1 1 

associated with the wave are u,, “ 0.2 m s and u., “ 0.05 ms, 

h V 

respectively. Similarly, the phase velocity Is 0.2 m s”K The vertical 
displacement caused by the wave, which follows from equations 7.1 and 
7.2, Is about 50 m at 44 km. In the Boussinesq approximation used 
impllcity here, the fractional variations In temperature 
obey 6T/T “ -6p/p [Spiegel and Veronls, I960]; substitution yields wave- 
associated temperature fluctuations of about 0.07 K at 44 km. 

Below 25 km, the scintillations at 13 cm during both immersion and 
emersion exceed the predictions of the gravity-wave theory. This 
observation and other considerations (see below) suggest that 
tropospheric processes below 25 km excite gravity waves which transport 
energy and momentum upward through the atmosphere. The vertical flux of 
energy due to gravity waves is given by 6p«u^, where 6p is the pressure 
field associated with the wave and the overbar denotes a wave-cycle 
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average [Llndzen and Holton, 1968], Using the approximations above, the 
upward flux of energy at 44 km Is roughly 2»10" W m" (2 erg cm" s"'^). 
The vertical flux of horizontal momentum from gravity waves Is given by 
PUj^'u^ [Gossard and Hooke, 1975, p.l08]. The net flux depends on a 
summation over the momentum components of each wave In the 

superposition, which Is not possible based on the present observations. 

-3 -2 -2 -2 

However, an upper limit of 3*10 N m (3 <10 dyne cm ) can be placed 
on the net upward flux of horizontal momentum by assuming that the 
contributions from the waves add constructively, which amounts to the 
assumption of an anisotropic gravity-wave source. 

Viscosity and thermal conductivity cause the dissipation of small 
scale waves, and Impose a lower limit on the permitted gravity-wave 
scale size. This minimum scale size increases at higher altitudes 
because of the increase In kinematic viscosity that accompanies a 
decrease In atmospheric density [Hines, 1960J. The effect of 
dissipation In Titan's atmosphere was evaluated using approximate 
formulas derived by Gossard and Hooke [1975, p. 218], the temperature- 
pressure profile in Figure 7.4, and measurements of the viscosity 
of Ng appropriate to Titan's atmosphere [Vargaftik, 1975, p. 458]. 
According to these calculations, gravity waves of length scale a few 
hundred meters to a few kilometers, which correspond to the size of the 
observed irregularities, can propagate essentially unattenuated to at 
least 150 km In Titan's atmosphere; this result provides a basis for 
discussing the importance of the gravity-wave transport of energy and 
momentum. 

This gravity-wave energy flux Is believed to be insignificant 
relative to the rate at which solar energy is absorbed in the 
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stratosphere or at the surface [private cornmunication, R. E. Samuelson, 

1982]. Hence, gravity waves will not affect the atmosphere's thermal 

structure below 150 km; only at considerably higher altitijdes would this 

energy flux become significant. In contrast, the absorption of momentum 

-3 -2 

above 150 km at a rate of 10 N m would accelerate the atmosphere by 
about 7 m per Titan day, which may be significant in sustaining the 
cyclostrophic zonal winds described by Hanel et al. [1981] and Flasar et 
al. [1981]. 

The net transport of momentum does not require anisotropic gravity- 
wave sources, and can occur under quite general source conditions when 
gravity waves propagate through a region of vertical wind shear. 

Briefly, waves propagating upward through a region of sufficient wind 
shear are either absorbed or reflected depending on whether their 
components of horizontal momentum are generally parallel or anti- 
parallel, respectively, to the mean flow. The theory of reflection and 
"critical layer" absorption has been discussed at length [Gossard and 
Hooke, 1975, p. 170; Bretherton, 1969a, b; Turner, 1973, p. 24-29 and 
37-39]; this same type of interaction, where gravity waves transport 
horizontal momentum from the troposphere, where they are generated, to 
the stratosphere, where the momentum is absorbed by the mean flow at a 
critical layer, has been proposed by Lindzen and Holton [1968] to 
explain features of the zonal flow in the Earth's tropical atmosphere 
[also see Holton and Lindzen, 1972; Wallace, 1^J73; Cadet and Tiet.elbaum, 
1979], Until more detailed observations are available, one can only 
speculate that gravity waves may supply momentum to Titan's zonal 
circulation through this type of Interaction, It should be notefi that 
some momentum transport process is required to maintain the excess 

angular momentum in Titan's atmosphere [Flasar et al., 1981]. 
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The preceding discussion preseiited evidence In support of a 
gravity-wave Interpretation for the refractive Irreglarltles In Titan's 
atmosphere above 25 km. Although other Interpretations are possible, 
each suffers from Inconsistencies with the observations. For 
completeness, and to provide Indirect support for the gravity-wave 
hypothesis, several alternative explanations and their limitations are 
summarized below. 

Minor constituents In Titan's stratosphere Include CgHg, CgHg, 

^3^8* ^3^4’ [Maguire et al., 1981; Kunde et al., 1981]. Figure 

2 of Maguire et al. [1981] indicates that the condensation level of each 
of these gases lies between about 55 and 85 km. As the upper boundary 
of condensation is expected to be abrupt, the condensation products 
could Introduce layers In refractivlty that contribute to the scattering 
observed In Titan's stratosphere. Alternatively, thin stable layers are 
observed commonly in the earth's troposphere and stratosphere under 
cloud-free conditions [Harper and Gordon, 1980; Gage and Balsley, 1980; 
Rottger, 1980], and could also contribute to scattering In Titan's 
stratosphere. However, In either case, stable stratification cannot 
account for the Inferred horizontal variations In refractivlty. 
Furthermore, the agreement between the gravity-wave theory and the 
observations in Figure 7.5 suggests that these scattering mechanisms are 
not important In Titan's upper atmosphere. 

Thin turbulent layers In clear air have been observed in the 
earth's lower stratosphere [R. K, Crane, 1980]. These turbulent regions 
develop when an atmospheric layer becomes unstable with respect to 
vertical displacement because of excessive wind shear, A related 
measure of atmospheric stability Is given by the gradient Richardson 
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number, Ri = gn/(du/bz)^ [Wallace, 1973]. Here, the numerator Is the 
square of the Valsala-Brunt frequency, the denominator Is the square of 
the wind shear, and RI > 0,25 Is generally regarded as the condition for 
stability [Turner, 1973, p. 100], In the Interval between 25 and 75 km, 
(gll)^^^ > 4*10"^ requiring oU/dz > 8 m s"^ km"^ for Instability. A 
wind shear of th-is size has not been confirmed by any observations -- 
turbulent layers are not expected to be present In Titan's atmosphere 
above 25 km. 

A hypothesis that the scattered signal originates from acoustic 
waves also fails to explain the observations. An acoustic-wave energy 
flux of about 1 W m"^ would be required to account for the parameters in 
Table 7.1 [for the appropriate theory, see Gossard and Hooke, 1975]. 
However, there Is no plausible energy source for waves of this 
magnitude. For example, a source in the boundary layer near the surface 
would necessitate remarkably efficient conversion of the incident solar 
flux of about 3 W m' [Samuel son, 1983], and would deplete the 
Insolation needed to maintain radiative equilibrium at the temperatures 
in the lower atmosphere [cf. Samuel son, 1983; Linda! et al., 1983; 
Eshleman et al., 1983]. Furthermore, Figure 7.6 shows two theoretical 
profiles cf <n^ >, a quantity that is proportional to as discussed 
above. One curve is appropriate to internal gravity waves, or buoyancy 
waves (see Section 1.3), and also appeared in Figure 7.5; the other 
curve corresponds to freely-propagating acoustic waves. The latter was 
calculated using formulas analogous to those in equations 7. 1-7. 3 
[Gossard and Hooke, 1975], and the measured stratification of Titan's 
atmosphere at Voyager 1 immersion (Figure 7.4). The buoyancy-wave 
profile agrees with the measurements above 25 km in Figure 7.5, but the 
acoustic-wave profile diverges significantly. 
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Figure 7.6 Two theoretical profiles for the variance of refractive 
Index. Each would result from a different type of atmospheric wave 
propagating through the vertical stratification In Titan's atmosphere at 
Voyager 1 Itranerslon. One curve Is appropriate to Internal gravity 
waves, or buoyancy waves, and also appeared In Figures 7.5a and 7.5b. 

The second curve corresponds to acoustic waves. The two curves have 
been scaled by constants to facilitate comparison. 


Lower Atmosphere (0~25 km). The temperature follows the dry 
adiabatic lapse rate for in the lowest few kilometers of the 
atmosphere [Lindal et al., 1983], Consequently, the weak scintillations 
from this region during immersion {see Figures 7.2 and 7.5a) might be 
the result of propagation through a cloud- free convective boundary 
layer. Detailed analysis of the signals scattered from this layer Is 
restricted by the short time duration and limited signal-to-nolse ratio 
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of the data. In addition, the strong refractivlty variations between 5 
and 25 km altitude tend to obscure the Irregularities between the 
surface and 5 km when observed through occultatlon techniques. More 
work Is needed on this problem. It Is Interesting to note that boundary 
layer turbulence and penetrative convection can excite gravity waves 
[Stull, 1976; Elnaudi et al., 1978; Nagpal, 1979; Harper and Gordon, 
1980] that would propagate Into the stratosphere. 

Between 5 and 25 km, the temperature lapse rate Is subadlabatic and 
apparently stable [Lindal et al., 1983; Eshleman et al., 1983], and the 
13-cm scintillations exceeded the predictions of the gravity-wave model 
(Figure 7.5). A comparison with the earth's atmosphere [Stull, 1976] 
suggests that penetrative convection Is unlikely to extend more than a 
few kilometers Into the stable region above the turbulent boundary 
layer. Consequently, the strong scattering In the middle troposphere 
perhaps was associated with propagation through stratified layers rather 
than dynamics. Lindal et al. [1983] and Eshleman et al. [1983] found no 
direct evidence for CH^ condensation in the radio occultatlon 
observations, and favor models In which the CH^ concentration Is at most 
a few percent. In this case, layers In Titan's troposphere similar to 
those observed In the earth's troposphere and stratosphere under clear 
air conditions [Green and Gage, 1980] could cause the strong radio 
scintillations in the absense of any clouds. However, Samuel son 
et al. [1981] proposed that CH^ clouds may provide the dominant opacity 
In Titan's upper troposphere, although Hunt et al. [1983] currently 
dispute this claim. On the basis of terrestrial cloud models [Mason, 
1971], the refractivlty variations resulting from layered clouds could 
explain some of the strong scintillations In this region provided that 
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sufficient material has condensed. At present, we cannot distinguish 
between these two Interpretations. 

Analysis of the strong Intensity scintillations (see Figures 7.2 
and 7,5) may provide additional information on the refractivlty 
structures In Titan's troposphere. In particular, It may be possible to 
Infer the morphology of the irregularities, as was the case In Titan's 
stratosphere, and thereby distinguish stable stratification from dynamic 
processes. Any information concerning the variations in refractivlty 
between 5 and 25 km altitude would lead to a better understanding of the 
source for stratospheric gravity waves, and the lowest few kilometers of 
the atmosphere as well. 
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CHAPTER 8 

SUMMARY OF RESULTS AND SUGGESTIONS FOR FUTURE WORK 


8.1 Summary of Results 

A new procedure has been developed for systematically analyzing the 
amplitude and phase scintillations observed during radio occultatlons by 
planetary atmospheres (Chapter 3). First, conventional methods of 
spectral analysis are applied In order to reduce the data and to model 
the statistics of the measurements. Then, the resultant amplitude and 
phase spectra are Interpreted through objective comparisons with the 
predictions of the well-developed weak-scattering theory. This 
technique yields least-squares solutions for parameters that 
characterize the underlying atmospheric scattering structures. The 
covariance matrix of the solutions Is calculated by propagation of 
errors, thereby providing the measurement uncertainties. This procedure 
was applied to data obtained during the dual-wavelength (3.6 and 13 cm) 
radio occultatlons of Voyager by Jupiter, Saturn, and Titan. 

For Jupiter and Saturn, only data pertaining to the upper 
Ionospheres were considered. For the theoretical calculations, the 
occulting Ionosphere was modeled as a thin screen containing a power-law 
spectrum of anisotropic scattering structures. Comparisons between 
measured spectra and the theory produced solutions for four parameters 
In the model for the Ionospheric Irregularities: (1) the axial ratio, 

(2) the angular orientation, (3) the power-law exponent of the spatial 
spectripa, and (4) the magnitude of the spatial variations In electron 
density. Scattering was observed from Irregularities of approximate 
size 2 to 140 km at Jupiter and 2 to 60 km at Saturn; In both cases, 
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the lower and upper limits result from noise limitations and the method 
of analysis, respectively, and do not represent either Inner or outer 
scales. Within this size Interval, the three-dimensional spatial 
spectrum obeys an Inverse power law with exponent varying from 3.0 to 
3.7; moreover, and the fractional variations In electron density range 
from 1 to 20%. All Irregularities appear to be anisotropic; the axial 
ratios, which vary with altitude and latitude, are as small as 2:1 at 
low altitudes In the equatorial Ionosphere, and as large as 12:1 at high 
latitudes. Magnetic field orientations in the Ionospheres of Jupiter 
and Saturn were Inferred from the direction of alignment of the 
anisotropic Irregularities. The occultatlon measurements conflict with 
the predictions of current Jovian magnetic field models, but generally 
corroborate models of Saturn's field, and may prove useful In deriving 
improved field representations. These results were discussed at length 
In Chapters 4, 5, and 6. 

During the occultatlon of Voyager 1 by Titan, Intensity and 
frequency fluctuations occurred on time scales from about 0.1 to 1.0 s 
at both wavelengths whenever the radio path passed within 90 km of the 
surface. This observation indicates the presence of variations in 
atmospheric refractlvlty on length scales from a few hundred meters to a 
few kilometers. As the Intensity fluctuations showed similar variations 
with altitude during Immersion and emersion, the small-scale atmospheric 
structure does not appear to be a localized phenomenon. The theoretical 
Intensity spectra show best agreement with the measurements when 
horizontal variations in refractivity are Included in the models for the 
scattering irregularities -- evidence for a dynamic atmosphere. Above 
25 km, the altitude profile of intensity scintillation closely agrees 
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with the predictions of a theoretical model In which Internal gravity 
waves propagate with no attenuation through the vertical stratification 
In Titan's atmosphere. These results substantiate the hypothesis that 
gravity waves underlie the variations of refractlvlty In Titan's upper 
atmosphere. A cloud- free convective region In the boundary layer near 
the surface could act as a source for the waves. The net upward flux of 
energy and an upper bound on the net upward flux of horizontal momentum 
due to gravity waves are approximately 2«10" W m" (2 erg cm’ s“ ) 
and 3-10 ‘ N m (3 *10 dyne cm ), respectively. Although this energy 
flux is not significant In determining the thermal structure of the 
atmosphere below 150 km, the vertical transport of momentum by gravity 
waves could be Important in maintaining the zonal stratospheric 
circulation. These results were described In detail in Chapter 7. 

8.2 Suggestions for Future liork 

8.2.1 Development of a New Technique for Analyzing Strong 
Scintillations. The theory derived in Appendix A is valid only for the 
case of weak scattering. However, periods of strong intensity 
scintillation were observed during the Voyager radio occultatlons by 
Jupiter, Saturn, and Titan (for example, see Chapter 7); a more general 
theory is required in order to interpret these data. In this regard, 
the development of new methods of analysis could improve the occultation 
technique by extending it into the strong- scattering signal regime. 

In a recent advance, a new theory was derived to predict the 
temporal spectrum of strong Intensity fluctuations during an occultation 
by an atmosphere with a finite scale height [Haugstad, 1982]. All 
previous derivations for the case of strong scattering neglected the 
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effect of differential refraction (or Fresnel zone distortion) 
associated with occultatlon measurements [cf. Rumsey, 1975], and 
therefore were Inappropriate for interpreting the strong Intensity 
scintillations mentioned above. A difficult numerical Integration Is 
required to Implement the new formula; after this has been accomplished. 
It will be possible for the first time to analyze strong Intensity 
scintillations during planetary occuHatlons. The theoretical spectra 
can be compared to the measurements by following the least-squares 
procedure described In Section 3.4. At present, one area of uncertainty 
In the potential of this new technique Involves the sensitivity of the 
strong-scattering spectra to the underlying dynamic structure In the 
occulting atmosphere. 

8.2.2 Further Analysis of Voyager Radio Occultatlon Data. Our 
knowledge of the neutral atmospheres of Jupiter and Saturn might benefit 
from analysis of the signal fluctuations observed during the Voyager 
radio occultatlons. Weak scintillations can be analyzed by following 
the procedures which were applied In Chapter 7 to Interpret data from 
the Voyager 1 occultatlon by Titan. Strong Intensity scintillations can 
be attacked by using the new technique described above. In this way, 
the structure and dynamics of these atmospheres can be Investigated on 
length scales near Fresnel-zone size (see Tables 4.1 and 5.1) and at 
pressures from about 1 mbar to more than 1 bar. 

Strong Intensity scintillations were observed in connection with 
propagation through Titan's lower atmosphere (5-25 km altitude; see 
Figures 7,2 and 7.5). Present uncertainties regarding the nature of the 
refractive Irregularities In this region might be resolved through 
application of the technique described above. By comparing calculated 
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theoretical spectra with the measurements, spatial characteristics of 
the atmospheric structures, such as the axial ratio and size 
distribution, might be inferred. (See Chapters 4, 5, and 7 for results 
derived from observations of weak scattering.) This characterization 
could lead to a better understanding of the structure and dynamics of 
Titan's troposphere, and might allow an identification of the source of 
the stratospheric gravity waves. 

While internal gravity waves in the earth's atmosphere have been 
studied extensively, the role of gravity waves in other planetary 
atmospheres remains uncertain. In fact, the Voyager occultation 
measurements at Titan constitute the only detailed observations of 
progressive gravity waves in the atmosphere of another planet. In 
previous stellar and radio occultations , the unambiguous identification 
of gravity waves as scattering sources has not been possible. For 
example. Woo et al. [1980a] Initially speculated that scintillations 
observed during the radio occultations of Pioneer by Venus were the 
result of trapped gravity waves in a vertically stable region of the 
atmosphere. However, Woo et al. [1982] subsequently revised this 
interpretation and attributed the scattering to turbulent layers 
generated by wind-shear instability. It also has been suggested that 
gravity waves cause the intensity fluctuations observed during stellar 
occultations by the upper atmospheres of Jupiter [French and Giorasch, 
1974] and Uranus [French et al., 1983]. Although the discovery of 
internal gravity waves at Titan lends credibility to this 
interpretation, the hypothesis remains speculative, since the underlying 
cause of intensity fluctuations in stellar occultations is at present 
the subject of controversy (for example, see the review by Elliot 
[1979]). 
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A study of atmospheric dynamics at Jupiter, Saturn, and Titan, as 
outlined .^bove, could promote a more general awareness of the role of 
nravlty-wave phenomena In planetary atmospheres. Progress In 
understanding the dynamics of these atmospheres could provide a basis 
for comparisons with the earth and with other planets as studied by 
occultation methods, thereby contributing to atmospheric studies at 
several planet?, 

8.2.3 Magnetic Field. Studies at Jupiter and Satu^ ’n. In a 
spherical harmonic magnetic field representation, the quadrupole and 
octupole terms In the expansion are relatively more Important In the 
ionosphere than at the remote position of spacecraft magnetometer 
measurements (see Section 6.5). Consequently, the field orientations 
Inferred from radio occultatlons could provide useful constraints on the 
higher-order coefficients In field models. For this reason, our 
knowledge of the internal magnetic fields of Jupiter and Saturn could be 
Improved through the development of methods for deriving field models 
from the combined magnetometer and radio occultation measurements. 

Two fundamental areas of research might benefit from future 
measurements pertaining to the magnetic fields of Jupiter and Saturn. 
First, Improved spherical harmonic field representations or evidence for 
secular variation could provide constraints for models of the Jovian and 
Saturnian Interiors and dynamos [Hide and Mai in, 1979; Stevenson, 

1980]. Second, efforts toward obtaining a more detailed field model 
near these planets might lead to a better understanding of the source of 
low-frequency radio emissions from both Jupiter [Warwick et al., 1979a, 
b; Boischot et al., 1981] and Saturn [Kaiser et al., 1980, 1981; Warwick 
et al., 1981, 1982], and an explanation for the modulation of the 
radiation at the planets' rotation rates. 

no 
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Observations of Faraday rotation In future radio occultatlons could 
provide the strength of the field component along the ray path 
[Eshleman, 1973]. A spacecraft orbiting Jupiter or Saturn and equipped 
for radio occultatlon measurements of both scintillations and Faraday 
rotation would be well suited for further study of the planetary 
magnetic field. When conblned with the conventional magnetometer data, 
a sequence of radio occultatlon measurements at a range of latitudes and 
longitudes could produce a more accurate field model near the planet. 
These types of measurements are planned as part of a future 
Investigation of the Jovian system using the Galileo ortiber [private 
communication, H. T. Howard, 1903]. 
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DERIVATION OF THEORETICAL FORMULAS FOR WEAK SCATTERING 
A.l Introduction 

The theory for propagation of electromagnetic waves through a 
weakly irregular refracting medium originated in the studies of radio 
star and satellite scintillations caused by the earth’s ionosphere [for 
example, see Booker et al., 1950; Newish, 1951; Fejer, 1953; Bramley, 
1954; and Ratcliffe, 1956]. The application of the weak-scattering 
theory to an occultation geometry has been discussed by Haugstad 
[1978b], Hubbard et al. [1978], and Woo et al. [1980b], among others. 
Although the results are not original, expressions for the theoretical 
covariance functions and scintillation spectra are derived here for 
completeness, and to provide a basis for the discussions in Appendix B 
(a limit for validity of the weak-scattering theory) and Appendix C {a 
physical interpretation of the weak-scattering theory). 

In Sections A. 2 and A. 3, the theory is derived for scattering from 
a phase screen with a spatially homogeneous average refractive index. 
Section A. 4 gives some convenient definitions for l'.e transverse 
coordinates in the phase screen. In Sections A. 5 and A. 6, the 
theoretical formulas are modified to account for the effect of Fresnel 
zone distortion (see Appendix D). Finally, Section A. 7 discusses the 
relationships among amplitude, log amplitude, and intensity 
fluctuations. 
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A. 2 Covariance Functions: Homogeneous Background 

A. 2.1 Approximate Forumlas for the Scattered Field. At the 
magnetic field magnitudes In Jupiter's ionosphere Inferred from 
spacecraft magnetometer measurements [Connerney, 1981] j the electron 
gyrof requency , f^, Is 10 to 20 MHz, while at the field strengths In 
Saturn's Ionosphere CConnerney et al., 1982], f “ 1 to 2 MHz. The 
observed Ionospheric Irregularities appear to be anisotropic and field 
aligned (Chapters 4 and 5), which suggests that the electron collision 
frequency v < f . In addition, the plasma frequencies, f , are less 
than 4 MHz and 1 MHz at the maximum electron densities measured In 
Jupiter's upper Ionosphere [Eshleman et al., 1979b, c] and Saturn's 
upper Ionosphere [Tyler et al., 1981a, 1982], respectively. 

The radio frequencies used for the Voyager occultatlon experiments 
(2295 and 8415 MHz) greatly exceeded f and f at both Jupiter and 

C p 

Saturn. For this reason, the effect of the planetary magnetic field on 
the refractive Index In the Ionospheres was extremely small [cf. Jordan 
and Balmain, 1968, pp. 687-698]. Nevertheless, a linearly polarized 
microwave signal propagating through Jupiter's ionosphere would be 
expected to experience appreciable Faraday rotation. However, the 
observation of this phenomenon was not possible with the Voyager 
spacecraft, which transmitted signals In right-hand circular 
polarization. 

In the remainder of this appendix, the phenomenon of Faraday 
rotation is neglected; the magnetic field then has no significant effect 
on direct propagation through the planetary ionospheres. (Note, 
however, that the magnetic field affects the morphology of the plasma 
irregularities, and thereby Indirectly Influences the scattered 
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signals.) The positions of the transmitter and receiver enter 
symmetrically into the expression for the scattered field (equation Al, 
below). Using this property, the roles of the spacecraft transmitter 
and earth-based receiver are interchanged in what follows, which 
simplifies the derivation and discussion without affecting the 
results. In particular, the electromagnetic wave impinging on the 
occulting planetary atmosphere can be treated as planar, since the 
distance from earth to the planet is greater by a factor of at least 700 
than the distance from planet to spacecraft for the occultations 
considered here. 

Figure Al shows the idealized experimental geometry. A 
monochromatic plane wave of unit amplitude propagating in the x- 
direction illuminates the occulting atmosphere, modeled as a screen of 
thickness L with transverse coordinates y, z. The spacecraft radio 
receiver is located a distance D beyond the midpoint of the screen, 
where D >> L. The refractive index in the scattering region is 
expressed as n(r) ~ 1 ■*" where n^ corresponds to the average 

electron density, and n^ represents the irregularities (<n^> = 0, where 
<> denotes a statistical average). For the ionospheric occultations 
considered here, the gradient in had no appreciable effect on the 
scattering; hence, in the model for tiie ionosphere, n^ is considered to 
be homogeneous { ^ 0 ^= 0 ). Furthermore, 1 + n^* 1 for the electron 
densities in the ionospheres of Jupiter and Saturn (see equation 2,1 and 
Figures 4.2, 4.5, 4.7, 4.9, 5.2, and 5.4). It is also assumed that the 
refractive variations obey the Taylor hypothesis so that they can be 
considered "frozen" on time scales of interest. During the occultations 
of Voyager, the spacecraft velocity was sufficiently large to ensure 
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that the only significant time variation resulted from motion of the 
raypath across an essentially fixed field of spatial variations in 
refractive Index (see Tables 4.1 and 5.1). 


PLANE 

WAVE 



Figure A1 Idealized experimental geometry and coordinate system. 

A monochromatic plane wave propagating in the x-direction illuminates 
the scattering volume modeled as a screen of thickness L with transverse 
coordinates y, z. The origin of coordinates is at the center of the 
occulting planet. A radio receiver is situated a distance D beyond the 
screen, where D >> L. 
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In the Born approximation, the electric field at the receiver Is 
tiic of the undisturbed Incident wave, and a small 

perturbation due to scattering, which can be expressed as [cf. 

Tatarskii, 1971, p. 108] 




///- 

V 




nt(r' ) ?g(r' ) 


d^r' 


» 


where 



(Al) 


Here, k Is the radio wavenumber, z Is the coordinate unit vector, the 
integration extends throughout the volume of the screen, V, and time 
dependence exp[-j2nf^t] has been suppressed, where f^ is the radio 
frequency. Approximations used in deriving equation Al Impose the 
following restrictions on Its validity [cf. Tatarskii, 1971, p. 103-106] 


k |r-r' I >> 1 and << X , 


(A2) 


where \ Is the free-space radio wavelength and Ji. Is the smallest 
irregularity which contributes significantly to the scattered field (as 
opposed to the smallest irregularity present). The above inequalities 
are easily satisfied for the Voyager occultatlons by Jupiter, Saturn, 
and Titan, where kjr-r' | > 10® and \/A < 10“®. [Values of Ji are 2 km 
for Jupiter (Chapter 4), 2 km for Saturn (Chapter 5), and 0.2 km for 
Titan (Chapter 7), while \ “ 10"^ km]. 


117 



Derivation of Theoretical Formulas 


In the present case, the scattering angles are of order \/A or 
smaller [Tatarskii, 1971, p. 115], Consequently, the principal 
contributions to arise In the region of the screen contained within a 
radius M)/A from the straight-line path connecting spacecraft and 
earth. For mathematical convenience, the limits of Integration over y 
and z can be extended Indefinitely {-“ < y,z < +®), since the remote 
regions of the screen produce a negligible contribution to f . 
Furthermore, the small-angle approximation allows the expansion 


If.f. I . (x-x. ) . . (A3) 

* 8 (x-x')^ 


Using A3, and neglecting second order corrections 


r,(?) 


^0^"^) -JT /JJ 


(A4) 


•n^(r' ) d^r' , 


Terms through second order in the expansion A3 have been retained in the 
exponential factor, which requires 

L(y-y')^+ « 2n . (A5) 

8(x-x* 

For the region of the screen where the scattered signal originates, this 
Inequality can be expressed as 
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where the left-hand side Is the radius of the Fresnel zone. For the 
occultatlons considered here, the condition Is easily satisfied (see 
Chapters 4, 5, and 7 for parameter values; typically, 5 km, 

A * 1 km, and k " 10 km). 

As the scattered field remains polarized In the same sense as the 
Incident field, both will be treated subsequently as scalar 
quantities. By convention, the weak-scattering theory solves for the 
first-order contributions to the log amplitude, Xi and phase, <p: 

E(r) = £<,(?) exp{x + j4») « E^(r) [l + x + jl*] . (A7) 

From expression A4, 


E'(r) H 


x(r) + j<)(r) 


k*" 

2n 


/// 

V 




L(y-y‘)^+ iz-z'f]} 


(A8) 




n^(r‘ ) d^r* 


A. 2. 2 Statistics of the Scattered Field. In a statistical 
description of the scattered field, the quantities of Interest are the 
autocovariance functions of amplitude and phase 

B^(ri,r 2 ) = <x(ri) x(r 2 )> 

(A9) 

= |Re L<E'(ri) E‘(r 2 )> + <E'(ri) E'(r 2 )>J . 
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B|jj{ri,r2) = «>{ri) 4>(r2)> 


= |-Re [<E'(ri) E'(r 2 )> - <E‘(ri) E'(r 2 )>j 


(AID) 


<ind the cross covariance of amplitude and phase at the same radio 
wavelength 

Bj^^(ri,r2) = <x(ri) 4>(r2)> 

(AU) 

= i Im [<E'(ri) E‘(r2)> ‘ <E'(ri) E'(r2)>] • 


In Equations A9-A11, Re and Im signify the real and Imaginary parts, and 
overbars indicate complex conjugates. 

It Is assumed that the spatial variations of refractive index In 
the phase screen can be described as a locally homogeneous random 
field. In this case, the spatial covariance of refractive Index, B^, 
depends only on the separation of the observation points and obeys the 
following Fourier transform relationship [Tatarskii, 1971; Yaglom, 

1962]: 


8^(C,ti,C) h <n^.(x+^,y+T],z+C) n^(x,y,z)> 

+<n 

= /// expLj(ic^e + KyTi + ®(x) dic^dx^dif^ , (A12) 

^G0 

$ (k) = — L_ /// expi-j(ic £ + X T) + X c)J B (^,Ti,C) d£ dn dC . 

” (2ii)‘^ X y z n 

Here, Is the three-dimensional spatial spectrum of variations in 
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refractive index, and k , k , and k are the Cartesian, components of the 

A jf 

■f 

wave vector, tc. 

The general expressions for the covariance functions in equations 
A9-A11 will be evaluated for the present case by using equation A8; two 
moments of the scattered electric field must be calculated. First, 

<E' (ri) E’ (ra)> » /// /// ) n^(r* ' )> 

Alt \t ]J ^ ^ 


L(yi-y')^+ (zi-z')^]} 


(A13) 


expl ^ -^ ’^^ -rr y [(y2-y‘ ' (Z2-z")'^]} 

. (X 2 ~ X ” ) 


,1 1 \2. 


,> < \2 I 


d^r'd^r' 


After substitution from A12, the integrals over y', z', y' ' , and 2 '' 
separate into factors such as {for y') 


+® 


1 2 /^exptjy' + ^y' > 


..^2 


(A14) 


which can be Integrated analytically to obtain 


Tl(Xl-X ) 1 yo X (XX"X' ) 

= (1 + j) ( j^'^exptjiiCyyi- J1 . 


W 


(A15) 


121 



Derivation of Theoretical Formulas 


After Integration over x' , x' ' , y' , y' ' , z' , and z' ' , expression A13 
reduces to 


where 


- Sln[(p + K ) L/2] , 

<E'(ri) E'(r 2 )> «= /// VV )/2 


•expijL-P(Xi-X2) + <y{yi- y2) + 22 )]} 

2 2 


(A16) 


Consider the two factors In the Integrand of A16 which depend 

on K , In the present study, It Is assumed that * obeys a power law of 
A n 

the general form (see section 2,3.1) 


$ tie •< 

n 




(A17a) 


where 2n/ic^ corresponds to the outer scale. On the basis of the results 

2 2 2 2 

In Chapters 4, 5, and 7, It was concluded that k + < + k >> k 

A Jr fc U 

for the Voyager observations at Jupiter, Saturn, and Titan, Then, to a 
good approximation 


^n 





<2)-P/2 
z ^ 


(A17b) 


2 2 1/2 

As a function of k , * has width of order (k + tc ) ' , and Is 
X n y z 

symmetrical about = 0, The second factor 


sinL(P + K ) L/2J j, 

^ I — K tT 72— ) 


(A18) 


122 



Derivation of Theoretical Formulas 


has area 2nL and width 2n/L. As 2ti/L Is extremely narrow r.ompared to 
the width of (the observed Irregularities are much smaller than the 
thickness of the phase screen), J behaves much like a delta function, so 
that the Integration of A16 over tc^ yields 




- 2irk^L // «^(0,Ky,ic^) 


(A19) 


•exp{j[-p(xi- X 2 ) + Ky(yi- ya) + ^ 2 )]] ' 


where approximation Is easily 
satisfied, which can be seen by considering A17b and the definition 
of P: 


2 2 2 
B + K+K 

y z 


(«y + V)U + 


2 2 
% *^z 

4k^ 


< t<y^+ 


(A20) 


t9 



9 


by condition A2. As Indicated by the condition = 0 in the argument 
of the field scattered at small angles is insensitive to refractive 
variations along the line of sight. 
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An expression for <E'(ri) E'(r 2 )> can be derived in a similar 
fashion by using equations AS and A12, and evaluating the integrals over 
X*, x'*, y', y", z‘, and z“ {as in A14 and A15): 

<E'(ri) E'(r 2 )V => -k^ /// «,(x) 

m 00 


•expljL-p(xi+ X 2 ) + <y(yi- yz) + K^(zi- Z 2 )J} 


s1nL(V S) L/2j sinL(K^- p) L/2j 

(V 


(A21) 


dK^dXydK^ 


Consider the x^-dependence of the integrand in A21, which contains 
a factor similar to J in A18 


_ sin[.(Kj^+ p) L/2j sin[(Kj^- p) L/2j 


(A22) 


Using the same arguments as applied previously to the function J, K can 
be approximated as a delta function when integrating over 


<E'{ri) E'(r 2 )> “ // ^^{O.Ky.x^) 


(A23) 


*exptji-p(xi+ X 2 ) + «y(yi- yz) + Zz}jl dfydx^ . 
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When substituted Into equations A9-A11, A19 and A2;3 yield 
expressions for the autocovariance and cross covariance functions of 
amplitude and phase: 


Bj^(ri,r2) 

B^(ri,r2) 


// «^(0 


.Ky.ic^) exptJi(c^{yi- yz) + zg)]} 


(A24) 

•lcos[(3(xi - X 2 )J + cosL^(x< + X 2 )]| dK^dK^ , 


Sj^^(ri,r2) “ If *n(0*''y*''z) exp(j|.iCy(yi- y2) + Z 2 )}} 

(A25) 

•ls1n[p(xi - X 2 )J + — s1n[p(xi + X 2 )J 1 dtc^dic^ , 

Note that *p(0ify so that the above 
expressions are real valued. Equations A24 and A25 are Identical to the 
results obtained for scattering from a slab of finite thickness In the 
Rytov approximation [Yeh and Liu, 1982; also see related comments 
concerning the Rytov and Born approximations In Appendix B]. 

Consider the parameter 3L In equations A24 and A25, By definition 

2 2 

If the Irregularities are larger than the radius of the Fresnel zone, 

1/2 

(XL) , as seen by an observer positioned near x = L/2, then 

pL << n and (sin PL)/pL « 1. Linder this condition, the field emerging 

from the scattering volume can be calculated accurately using 
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geometrical opticsi and tha occulting atmosphere behaves like a thin 
phase-changing screen CSal peter, 1967; Bramley, 1977]. 

Noting that x^- Xg < L for the observations reported here, and 
applying the thin-screen approximation 


Sj^(n,C) 




Jj expLj(« Ti k c)J 






I 


(A27) 




B^^(n.C) » 2nk^L // ®n(0.iCy .ic^) expLj(K^Ti + ic^c)j 


2 2 

•s1n[-|jpj cos[-|j~J « 


(A28) 


2 2 2 

where k ° and Xj+ Xg * 2D. For a homogeneous spectrum, 

the covariance functions depend strongly on the transverse separation of 
the observation points n = y^- and Q = z^- Zg, while the weak 
dependence on the separation along the line of sight, x^- X 2 > can be 
neglected. 

Let us assume that the amplitude of the Irregularities Is a 

constant fraction of the background refractive Index. The quantity 

(2rH)^^^ Is then a good approximation to the phase screen thickness L, 

where r Is the planet's radius and H Is the atmosphere's scale height. 

4 3 

For Jupiter, r « 7»10 km, and H « 10 km in the upper ionosphere 

[Eshleman et a1., 1979b, c], so that L « 10^ km; for Saturn, 

4 3 

r » 6“10 itm, and H « 10 km in the upper Ionosphere [Tyler et al., 
1981a, 1982], so that L » 10^ km; In Titan's atmosphere near the 
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tropopause, r • 2600 km and H ^ 16 km [Tyler et a1. . lysia], which give 
L " 300 km. At the wavelengths of the Voyager data and the values 
of A In Sections 4.3, 5.3, and 7.3, pL < n/4 so that {sin pL)/pL - 1. 
Therefore, the thln-screer approximation can be applied safely when 
considering the signal scintillations measured during the Voyager 
occultatlons. 


A. 3 Temporal Spectra; Homogeneous Background 

The temporal covariance of amplitude fluctuations, R (tj. Is 
related to the spatial covariance function 1n A27 through 




(A29) 


where v^ and v^ are the transverse components of spacecraft velocity 
(i.e., parallel to the phase screen), and t Is the time delay of the 
measurement. By definition, the temporal Fourier transform of A29 
yields the spectrum of amplitude fluctuations 

1 

w^(u) = J Bj^(v^t,v^t) exp(-jo>T) dT , (A30) 


where w Is the fluctuation frequency. After substitution from A27 and 
integration over x (which produces a delta function) 


+“ 


w^(u)) = 2nk L // * (0,K ,K ) 6 Lk V + K V,- uj 


n' ’ y* z' *■ y y z z 


.2, A 


•sin L^J dic^dK^ 


(A31) 
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For the case of a homogeneous average refractive index, the transverse 
raypath velocity in the phase screen, v*, is given by (see Appendix 0) 

^ A A 

V* “ ''yV + v^z . (A32) 

From this consideration, the physical significance of the delta function 
in A31 becomes clear: under the Taylor hypothesis, the contributions to 

amplitude fluctuations at frequency w arise from irregularities with 
wave vectors that satisfy 


w = , (A33) 

where 

V V * V • 


Changing variables of Integration, A31 becomes 


2 ” 

WyC***) " 2itk\ / J (0,KCOse,Ksine) 6[Kv.*cos(e-Y) - «J 
^ 0 0 ^ 

(A34) 

•sin^[-!^J K dK de . 


where k^= Kcose, KSine, and y is the angle between v* and the y- 

axis. Figure A2 shows the coordinate system used in equation A34 along 

with the angles e and y. | 
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Z 



Figure A2 Phase-screen coordinate system. The angles 
Y and e indicate the orientations of the transverse raypath 
velocity, v*, and the transverse wave vector, icj^, respectively. 
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Equation A34 can be integrated by using the following property of 
delta functions 


+® 


/ 6Lf(x)j dx = 



(A35) 


where f(x^) ° 0 and f (x^) ^ 0. lioting that the argument of the delta 
function in A34 has value zero only if k > w/v*. and applying A35 to the 
integration over e 


w (u) = 2i*^L / - 
^ w/v* ( 


TT72 ^ «.(0,KCOSe. .Ksine. 

2 

•sin^L-^J X dK , 


) 


(A36) 


^1’ 


Y + COS 


‘^[(i)/KV*j 


EgS Y - COS"^[(i)/lCV*) 


The lower limit for ic- integration is a consequence of condition A33. 

Similar expressions can be derived for the phase and cross 
amplitude-phase spectra 


2 +“ . 2 

w (w) = 2Tik L / - "2 ' '2 2TT72 ^ ^_(0,KCOse. ,<sine. ) 

^ (D/v* (ic^v/- i=l " ^ ^ 


2 1 I . 

“COS L~2jrJ X die , 


(A37) 


w .(w) 


= 2vk‘ 


+<» 2 

,22^ 2^172 il ^ 

ai/V* (k v/- u i = l " ^ ^ 


(A38) 


1 . ik'^D-i . 

•5- sinL-T;;~J X dK . 
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A. 4 Transverse Coordinates in the Phase Screen 

In order to complete the theoretical discussion of weak scattering 
from a phase screen with a homogeneous background, a specific definition 
of the transverse coordinates (y, z) is required. A simple choice for 
the orientation of the y- and z-axes follows directly from the 
fundamental characteristics of the refractive irregularities. In the 
ionospheric occultations considered here, the magnetic field influenced 
the shape of the plasma irregularities, and hence the characteristics of 
the scattered signals.. In contrast, the overall effect of the 
ionosphere's finite scale height was negligible. Consequently, in a 
convenient alignment of the transverse coordinates, the y-direction 
coincides with the projection of the magnetic field onto the y-z plane 
(see Section 2.3.1 and Figure 2.1). This coordinate system was adopted 
in Chapters 4, 5, and 6 when applying equations A36-A38 to interpret 
scintillations observed during the occultations of Voyager by Jupiter 
and Saturn. 

In Sections A. 5 and A. 6, equations A27, A28, and A36-A38 are 
modified to account for the effect of Fresnel zone distortion. The 
resulting formulas were applied in Chapter 7 to interpret weak 
scattering from Titan's neutral atmosphere during the radio occultation 
of Voyager 1. As a consequence of atmospheric stratification, the 
equations now are expressed most conveniently by using a coordinate 
system in which the y-axis coincides with the visible limb of the 
planet, and the z-axis is directed vertically, opposite to the 
gravitational acceleration at ray periapsis (see Section 2.3.2 and 
Figure 2.2). In general, this y-z orientation differs from the field- 
aligned system defined above. 
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A. 5 Covariance Functions: Inhomogeneous Background. 

In the preceding derivation, the average refractive Index In the 
phase screen was assumed to be essentially homogeneous. In order to 
remove this restriction from the previous results, the refractive index 
now Is expressed as n(r) = 1 + n^(r) + n^(r), where np(r) represents the 
spherically symmetrical quiescent atmosphere. The covariance functions, 
A27 and'A28, can be generalized accordingly by noting that the 
Inhomogeneous background affects the scattering In two ways {see 
Appendix D). First, for an observer at the spacecraft, the Fresnel zone 

Is distorted to horizontal and vertical dimensions {(f 

1 / 2 ’ ^ 
and ' , respectively. Second, horizontal and vertical 

separations 6y and 6z In the phase screen are transformed by 

differential refraction to equivalent separations and 6z/<j»^ in the 

plane of the spacecraft (x = D in Figure Dl). Expressions which account 

for the Inhomogeneous background are obtained by Introducing these 

corrections into A27 and A28: 
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B 


X4> 


+ CO 


(ti,C) = 2rt<'^L // 4 (0.< 

y 


“sinL^(0^ 


2 

. K + 4) K 

y y z z 


.^2) exp[j(<^y(CyTi + 


(A40) 
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Two Important properties can be Inferred from equation A39 for the 

case of scattering from a neutral planetary atmosphere. First, 

when obeys an inverse power law with exponent p (cf. equation AI.7), 

it can be shown through a change of variables that B (0,0) = <x > has 

radio wavelength dependence [Tatarskii, 1971; Young, 1976]. 

2 

Second, from equation A12, <n^ >, and It Is clear that 

2 2 

<x > * <n^ >. Furthermore, as refractivlty Is proportional to density 
for nonpolar gases [Essen and Froome, 1951; Tyler and Howard, 1969; 
Eshleman, 1973; also see Section 2.2], both <n^ > and <x > are 
proportional to the variance of fluctuations In gas density 1n the 
occulting atmosphere. 

A. 6 Temporal Spectra: Inhomogeneous Background 

Temporal spectra also can be derived from equations A39 and A40 by 
following the procedure In Section A. 3, which yields 


w 


2 +® 1 

(id) = 2uk L / p 7, I « (0,KCOse, ,Ks1ne. j 

«D/v* (tc'^v/- 1=1 " ' ’ 

2 

*s1n^i-^((t)yCos^e^+ ii>2Sln^e^}] k die , 

e^= Y + cos"^ (w/kv*) , 6^= Y - COs"^((d/kv*) 


(A41) 
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(A43) 


In this case, differential refraction modifies v*, the raypath velocity 
In the phase screen, according to (see Appendix D) 

V* = fyVyy + <f^v^z . (A44) 

Equations A41-A43 reduce to equations A36-A38 when it>y= 1, which 
corresponds to a homogeneous average refractive index, n^(r) = constant. 


A. 7 The Relationship Between Amplitude and Intensity Fluctuations 
This section discusses the simple relationships among amplitude, 
log amplitude, and intensity fluctuations for the case of weak 
scattering, and summarizes the notation for these quantities. 

In equation A7, the electric field at the receiver was expressed as 

E(r) EQ(r)oexp(x + j4>) “ Eo(r)*U + X + . (A7) 

where represents the undisturbed field. The approximation is 

valid under conditions of weak scattering, since the fluctuations in log 
amplitude, x, and phase, ip, then satisfy x» (p << 1 (see Tatarskii [1971, 
pp. 253-258] and Appendix B). 
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The first-order Intensity follows from equation A7: 

J(r) " lQ(r)’U + 2x] . (A45) 

Then, from equations 3.1 and 3.4, the normalized amplitude, a(t), and 
the normalized Intensity, 1(t), are related to x(t) through 


a(t) “ x(t) , 
i(t)=2x(t) . 


(A46) 


By applying this result, the covariance function for Intensity can 
be expressed as 

B^(ri,r 2 ) = <1(ri) i(r 2 )> = A*B^(ri,r 2 ) . (A47) 

In particular, 

<1^> = 4»<x^> . (A48) 

Similarly, the temporal spectrum of intensity fluctuations is given by 
(u)) = 4<>Wj^(a)) . (A49) 

This completes the theoretical derivations. 
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A LIMIT FOR VALIDITY OF THE WEAK -SCATTERING THEORY 

In the Born approximation, only first-order terms are retained In 
the perturbation expansion for the scattered field. Consequently, the 
results In Appendix A remain valid only for small signal fluctuations. 

In discussing the condition for validity of the weak-scattering theory, 
It Is useful to define the Intensity variance 

= < L ' • (Bl) 

2 

where I = E^ »exp(2x) is the signal Intensity (see equation A7), and < > 

represents an average over time. For weak scattering, this quantity is 

2 2 

related to the variance of log amplitude through « 4<x > (cf. 

Section A.7). 

2 

An upper limit on consistent with the Born approximation can be 

established from the predictions of theories that properly account for 

strong scattering. In the weak-scattering theory, an Increase In the 

variance of refractive Index (the parameter Q In Section 2.3.1) scales 

the Intensity spectrum but does not change its shape on a logarithmic 

display. Marians [1975] has computed Intensity spectra based on the 

theory for strong scattering from a thin screen [Rumsey, 1975]; as shown 

in Figure 6 of Marians [1975], the shape of the Intensity spectrum 

begins to depart from the predictions of the weak-scattering theory 
2 

when Increases beyond about 0.3, Similar computations by Rino 

[1980] are consistent with this limit (see his Figure 5). Apparently, 

2 

the weak-scattering theory Is valid only when < 0.3, for which 
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<X^> i 0.08. For all data analyzed In Chapters 4, 5, and 7, < 0.25, 

so the weak-scattering theory can be applied safely. 

An empirical test of the vveak-scatterlng assumption was obtained as 
follows. Spectra of the Intensity fluctuations can be computed by 
normalizing the data In two ways 




I - <1> 
<I> 


Ig = ln{I) - <ln(I)> 


(B2) 


For weak scattering, x Is small so that exp(2x) ■ 1 + 2x, and the 
intensity spectra obtained from 1j and 1g are essentially Identical. 

When the limits of weak scattering are exceeded, higher-order terms must 
be Included In the expansion of exp(2x), and the spectra computed 
from 1j and 1g are different. A comparison of Intensity spectra based 
on these two normalizations confirms that the Voyager data considered 
here correspond to weak scattering. 

It should be noted that some authors have adopted a less stringent 

condition for the validity of the weak-scattering theory. Tatarskii 

[1971, p. 258] stated that the Rytov approximation Is justified 
2 

when <x > < 0.64. However, at that time, the Rytov approximation was 
believed to constitute an Improvement over the conventional perturbation 
techniques based on the Born approximation, which Is strictly limited to 
single scattering. More recently, It has been recognized that the Rytov 
approximation does not properly account for multiple scattering and 
actually is equivalent to the single scattering formulations, and 
limited to about the same range of validity < 0.3 [Fante, 1975; 


Validity of the Weak-Scattering Theory 

Brown, 1966]. Woo and Yang [1976, 1978] adopted Tatarskii's criterion 

and applied the weak-scatterIng theory to Interpret amplitude spectra 
2 

when <x > was as large as 0.37, which is well above the limit for weak 
scattering Imposed here. 

The preceding discussion strictly applies only to the case of 

scattering from a thin phase screen with a homogeneous average 

refractive Index. However, it can be shown [cf. Haugstad, 1982] that 

2 

approximately the same upper limit on <x > Is appropriate when the phase 
screen includes an inhomogeneous background. 
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APPENDIX C 


A PHYSICAL INTERPRETATION OF THE WEAK -SCATTERING THEORY 
C.l Introduction 

The general behavior of weakly-scattered, short-wavelength radio 
waves has been discussed previously LRatclIffe, 1956; Lee and Harp, 

1969; Briggs, 1975; Uscinski, 1977]. Here, those results are extended 
In order to gain some Insight Into the mathematical formulas of Appendix 
A. In particular, this appendix explains the basic relationship between 
the morphology of the scattering Irregularities and the appearance of 
the temporal spectra of the received signal. It 1s emphasized that the 
spectra are sensitive to variations In both the axial ratio and the 
orientation of anisotropic Irregularities. 

Section C.2 provides a basis for subsequent discussions by 
analyzing a simple scattering problem. Section C.3 then applies the 
results to account for the principal features of the amplitude and phase 
spectra. Certain characteristics of the cross amplitude-phase sp'ictrum 
are discussed briefly In Section C.4. These results have been published 
elsewhere [Hinson and Tyler, 1982]. 


£a.ge uiiank naa: sulmed 
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Physical Interpretation 


C.2 The Field Scattered from a Homoaeneous Phase Screen 

It Is useful to begin by analyzing a simple scattering prc.Mem; 
Figure A1 shows the geometry. A monochromatic plane wave propagating In 
the x-directlon Illuminates the scattering volume, modeled as a thin 
phase screen with transverse coordinates y, z. At position x » D, a 
spacecraft receiver samples the perturbed wave. Consider thd case In 
which the refractive Index In the thin screen Is given by 


n(x*,y',z') “ 1 + OjCos{Kyy' + k-^z') , (Cl) 

where the constant n^ << 1, and the average refractive Index Is 

Independent of position. Note that 2 «/k and 2n/x represent the 

y ^ 

characteristic length scales of refractive variations in the y- and z- 
dlrectlons, respectively. For short radio wavelengths, the radio 
wavenumber k >> <^. Far beyond the screen, at location 

(x = D, y, z), the scattered field, E^, will consist of two plane waves 
Inclined at small angles relative to the original direction of 
propagation [Ratcliffe, 1956]; 


Eg(x,y.z) « jMiexpLj(tCyy + k^z)] + expL-j(y + 
•expijLx(k^-K^)^'^^- 2nf^tJi • 


(C2) 


Here, the real constant M Is proportional to n^ and represents the 
strength of the scattered field, f^ is the radio frequency, and 
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By expanding 




k 



(C3) 


and supressing the propagation factor, exp[j(kx-2itf|.t) J, the amplitude 
and phase perturbations (x and respectively, in the notation of 
equation A7) can be expressed as ' 

2 

X = 2M»cos(.tc^ + K 2 ^)*sin(||^) 

(CO 

2 

= 2M»cos(iCyy + 

The pair of scattered plane waves interfere with one another in such a 

way that at locations where the amplitude perturbation is 

maximum (^for example, x = ^ — ) , the phase pert Tbation is zero, and at 

2 

locations where x is zero (for example, x=0), the magnitude of 4> is 

maximum. The size of the irregularity, 2 h/k, relative to the radius of 

the Fresnel zone, (\x)^^^, determines the strength of the interference 

2 2 

effect through the factors sin(-^l^) and cos(-g^). These factors, which 

appeared in equations A27 and A28 for the covariance functions, will be 

referred to as the Fresnel filters in amplitude and phase, respectively 

[cf. Briggs, 1975]. Note that amplitude fluctuations develop only at a 

2 

distance greater than about x « ti k/2< beyo.nd the screen. Conversely, 
at a given distance beyond the screen, only irregularities on the order 
of Fresnel-zone size or smaller, such that ic > it/(\x) ' , will cause 
significant amplitude fluctuations (cf. equation C4}. As the variation 
with X is gradual (k » k), the time variations in x and <t> result 
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essentially from spacecraft motion parallel to the phase screen. In the 
case above, the fluctuation frequency Is id = <»v*, where v* Is the 
transverse raypath velocity and I< Is the wavevector with 
components and (cf. Section A. 3). 

In order to model a planetary atmosphere, the phase screen must 
Include a distribution of Irregularities, such as the one In equation 
Cl, with various wave vectors as specified by the spectrum, (see 

equation A12 and Section 2,3.1). However, for a statistically 
homogeneous irregularities with different wave vectors are 

Independent [Tatarskii, 1971, p. 22-25; Ysglom, 1962]. In this case, 
the weakly-scattered electric field (see equations A27 and A2S) consists 
of a linear superposition of Independent components, each In the form of 
a pair of plane waves, such as the one just analyzed. This 
representation of the field as a spectrum of plane waves is useful for 
interpreting the temporal spectra of signal fluctuations [cf. Uscinski, 
1977]. 

C.3 The Temporal Spectra of Amplitude and Phase 

The principal features of the temporal spectrum of amplitude 

fluctuations, w (u), can be explained In terms of the simple concepts 
X 

discussed above. First, consider scattering from a power-law spectrum 
of isotropic irregularities (e = f = 1 in equation 2.4), so that 

*^(ic) « k'P . (C5) 

For the observations by Voyager, 3 < p < 4 (see Chapters 4, 5, and 7). 
The amplitude fluctuations at frequency uj originate from those 
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refractive Irregularities where the wave vectors satisfy k'v*^ u. The 
Fresnel filter, plus which givf.s the strength of the 

Irregularities, determine the size of each contribution to the scattered 
signal. As rapidly decreases with Increasing k, the principal 

contribution to w (w) arises from the "first pass band" of the amplitude 

A 

Fresnel filter, which corresponds to Irregularities of size near the 
Fresnel zone radius {see Figure Cl). This is a well' known result of the 
wave-optics theory [Tatarskii, 1971], Fresnil-zone-slze variations in 
the screen contribute to w (u) at frequencies roughly less than or equal 

A 

to 


✓2nlv*| 



(C6) 


depending on their orientation with respect to v*; when p = 3.5, 
w (w) Is nearly constant for w< w . At frequencies u » u (smaller 

A V v» 

scattering structures), w (u)} generally decreases as when * (if) 

A ' ^ 

obeys a power law (see equation 8 of Hinson and Tyler [1982]). The 
break frequency, w , Is often referred to as the Fresnel knee. Figure 
C2 shows w (w) for scattering from Isotropic irregularities; on a log- 

A 

log plot, the high frequency asymptote is linear, with slope 1-p. 
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Figure Cl Functions that determine the efficiency of refractive 

irregularities for producing amplitude and phase fluctuations. The 

spectrum specifies the strength of the irregularities with wave 

vector ic; here, * « p = 3.5, and the irregularities are 

isotropic. The functions sin^ [tc'^D/2k ) and cos (k D/2k) represent the 

Fresnel filters for amplitude and phase, respectively. For a refractive 

structure with length scale 2 ti/k, the contribution to amplitude or phase 

fluctuations is proportional to the product of and the appropriate 

Fresnel filter (see equation A27). The amplitude fluctuations arise 

primarily from irregularities within the "first pass band" of the 

1/2 

Fresnel filter for amplitude, near k = /2it/(?^D) . However, the 

Fresnel filter for phase has value unity for small k. Hence, the 

principal fluctuations in signal phase are the result of scattering from 
the largest refractive structures present in the atmosphere. 
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Figure C2 Theoretical amplitude (upper panel) and phase (lower 

panel) spectra for weak scattering by a power-law phase screen; all 

scales are logarithmic. These spectra were computed numerically using 

equations A36, A37, and 2.4. In each panel the upper curve applies 

when ‘^p('f) Is isotropic, while the lower curve corresponds to scattering 

by anisotropic irregularities with an axial ratio of 20:1. The two 

curves In each panel have been shifted vertically by 8 dB to accommodate 

comparison. In the isotropic case, the Fresnel knee in w (w) occurs at 

frequency below which the spectrum is nearly constant. For the case 

of highly anisotropic irregularities, the Fresnel knee is reduced in 

frequency to and at lower frequencies the scintillations are 

suppressed relative to the isotropic case. Note that an increase in the 

axial ratio results in a shift in the frequency of the Fresnel fringes 

in both panels. For these calculations, jv,^j - 7.07 km s~^, 

Y = 30°, A = 0°, p = 3.5, and (XO)^^^ = 5 km. From equations C6 and 

C7, w = 2n and a>' = it, 
c c 
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Next consider scattering from highly anisotroic irregularities 
(e » f in equation 2.4). In this case, the phase screen has refractive 
variations in only one direction -- all wavevectors, k, in are 

essentially parallel. Consequently, temporal amplitude fluctuations at 
the receiver result only from the component of raypath motion parallel 
to ic. On the basis of the discussion above, a maximum in w (u) due to 

A 

Fresnel-zone-size structures occurs fsear frequency 


w: 


✓2nv; 


(C7) 


where v^ is the component of v* which is parallel to ic. For magnetic- 
field-aligned irregularities at angle y to the motion of the raypath, 
vi = |v*|*siny (cf. Section 2.3.1 and Figure 2.1). Refractive 
irregularities that contribute to amplitude scintillations at 
frequencies less than u»' have dimension greater that (\D) . Due to 

the Fresnel filter, structures of this size are not effective in 
producing amplitude fluctuations, so w (ui) is suppressed relative to the 
isotropic case for w < At frequencies greater than w^(u) 
again generally decreases in proportion to The anisotropic 

spectrum also appears in Figure C2, 

It is important to appreciate the difference bet‘/een the amplitude 
spectra in Figure C2. In the Isotropic case, the Fi asnel knee occurs at 
a frequency below which the spectrum is nearly constant (for p = 
3.5). For the case of highly anisotropic irregularities, the Fresnel 
knee is shifted in frequency to u', and at lower frequencies the 
scintillations are suppressed. The amplitude spectra therefore depend 
on the two-dimensional spatial structure of refractive variations in the 
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phase screen; this sensitivity can be explained physically. As 
demonstrated above, the Fresnel filter, which results from diffraction 
effects and experimental geometry, introduces a preferred length 
scale, for scattering structures. With a knowledge of the 

geometry, in particular v* and 0, and a measurement of the frequency of 
the Fresnel knee, w' one can infer y. the principal orientation of the 

C 

irregularities in the anisotropic case (see equation C7). Occultation 
observations along a flyby trajectory thereby yield two-dimensional 
infonnation concerning atmospheric spatial structure. Furthermore, as 
demonstrated by the curves in Figures C2 and C3, w (w) is sensitive to 

A» 

the degree of anisotropy. Therefore, accurate measurements of the axial 
ratio do not require simultaneous observations using a spaced array of 
radio receivers, contrary to what some authors have claimed [see e.g. , 
Woo and Yang, 19'*8], 
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0.27T 27T 20fr 

FREQUENCY, rad/sec 

Figure C3 Theoretical amplitude (upper panel) and phase (lower 
panel) spsctra for weaK scatter! nci by a power-law phase screen; all 
scales are logarithmic. These spectra were computed numerically using 
equations A36, A37, and 2.4. In each panel, the top curve applies 
when is isotropic, the middle curve represents scattering by 

slightly anisotrcp'fc irregularities with an axial ratio of 1.7:1, and 
the bottom curve corresponds to an axial ratio of 3:1. Again, the 
curves have been shifted vertically by a few dB to accommodate 
comparison. Notice that an increase in the axial ratio enhances the 
Fresnel fringes at <»> > u' and suppresses che amplitude scintillations 
for u < For these calculations, (v^i = 7.07 km s , y = 90°, 

A = 0°, p = 3.5, and = 5 km. In this cast*, = 2it 

since y = 90° (see equations C6 and C7). 
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Oscillations, or "Fresnel fringes", occur at frequencies <*)>(«)' In 
each of spectra In Figures C2 and C3. A comparison of curves In 
Figure C3 shows that the Fresnel fringes are enhanced by anisotropic 
Irregularities. Once again, this phenomenon can be explained on the 
basis of the arguments given above [cf. Briggs, 1975]. In the case of 
strong anisotropy, all wavevectors in are nearly parallel. 

Therefore, all Irregularities that satisfy w = <*v*, and contribute 
to w (u) at this frequency, have nearly the same magnitude and 
experience the same amount of Fresnel filtering. The amplitude Fresnel 
filter (see Figure Cl) strongly modulates w (u), producing local maxima 

A 

at frequencies 

w = (m)^^^u.’ , (C8) 

for m = 1, 3, i), ..., with alternating deep nulls (m = 2, 4, 6, ...). 

For the case of isotropic Irregularities, the situation is modified; 
there Is a range of wavevectors with different magnitudes and 
orientations contributing to w (u) at each w. These various wavevecv,ors 
experience different amounts of Fresnel filtering depending on their 
magnitudes. Consequently, in the summation of contributions at each a>, 
the effect of the Fresnel filter tends to be smoothed out; the resulting 
fringes are less pronounced than In the anisotropic case. 

Similar considerations apply to Interpretation of the theoretical 
phase spectrum, w^(u). As the Fresnel filter for phase has value unity 
for small < (Figure Cl), no Fresnel knee Is present In w^(i*>). The 
spectrum generally decreases as with Fresnel fringes similar to 

those In w (w), but at offset frequencies. Theoretical results 
for w^(u) also appear in Figures C2 jnd C3. 
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To summarize the discussion above, the theoretical model for 
refractive Irregularities In the Ionosphere contains five free 
parameters (Section 2.3.1); (1) the Inverse power-law exponent., p; 

(2) the axial ratio, e/f; (3) the orientation, y; (4) the strength, Q; 
and (5) the wave vector corresponding to the outer scale, As did 
not influence any of the observations In Chapters 4, 5, and 7 (see 
Section 3.4), this parameter was not considered here. (The theoretical 
curves in Figures C2, C3, and C4 were calculated In the limit 
where x^ -*■ 0; this parameter then had a negligible effect on the shape 
of the spectra over the frequency range of Interest). The remaining 
parameters affect the spectra as follows: the exponent, p, determines 
the slope (on a logarithmic plot) of the high-frequency portion 
both w (ii>) and w (oj). For anisotropic Irregularities, the 

K T 

orientation, y, strongly determines the frequency of the Fresnel knee 

In w (w), while the axial ratio, e/f, influences the shape of w (ui) at 

low frequencies. In addition, both e/f and y affect the detailed 

structure of the Fresnel fringes. Finally, the factor Q scales both 

w (oj) and w (<a) for weak scattering, but does not change their shape on 
X T 

a logarithmic display. 

As a direct consequence of these distinctive effects, the four 
principal parameters can be obtained from analysis of each measured 
amplitude or phase spectrum by following the procedure In Section 3.4. 
The strength of the determination depends on the values of the 
parameters themselves and the experimental conditions. Analysis of the 
solutions for p, y» and Q shows that for the Voyager observations, these 
parameters were well determined in a physically meaningful sense. The 
parameter e/f is best determined from measurements at low frequencies 
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(u) < (*>') which require long time segments of data. The relatively large 
uncertainties In e/f reported In Chapters 4, 5, and 7 are the result of 
the limited duration of the data from any one altitude region In the 
occulting planet's atmosphere. 

The preceding discussion can be generalized to the case where the 
refractive background Is Inhomogeneous; the effects of the free 
parameters In on the scintillation spectra are essentially 

unchanged. 

C.4 The Cross Amplitude-Phase Spectrum 

The cross spectrum, w^^(w), compares the magnitude and phase of 
correlated Fourier components In the amplitude and phase spectra of the 
received signal (see Sections A. 2 and A. 3). By the reasoning In Section 
C.2, at each frequency oi, consists of Independent contributions 

from irregularities with wave vectors that satisfy ic«v*= w. The 
behavior of each contribution can be understood by considering the 
result In equation C4 for scattering from the cosinusoidal phase 
screen. In that case, fluctuations In amplitude and phase caused by 
spacecraft motion parallel to the screen are either in phase or shifted 

by n radians, as determined by the signs of the Fresnel filters, 

2 2 

sin(K 0/2k) and cos(k D/2k). The relative size of the fluctuations also 
depends on these two factors. Clearly, since the phase of each 
contribution to Is 0 or it radians, w^^(<d) is real for scattering 

from a statistically homogeneous medium. 

Figure C4 shows the magnitude and phase of ;:s calculated 

from equations A38 and 2.4. Fluctuations at low frequencies result from 
irregularities with large dimension (small |k|) so that s1n(K D/2k) and 
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Physical Interpretation 


cos(tc^D/2k) are both positive and the cross spectnw phase Is 0. At a 
higher frequency, there Is an Increased contribution from length scales 
where s1n{K^D/2k) > 0, cos(ie^D/2k) < 0, and the phase of changes 

abruptly to n. In this way, the phase Is predicted to alternate between 
values of Ci and n, with shifts occurring at frequencies u = (m) ' for 
Integer values of m. Each shift In phase Is accompanied by a null In 
the magnitude of w^^(u)). 
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Physical Interpretation 




Figure C4 Theoretical magnitude (upper panel) and phase (lower 
panel) of the cross spectrum, Wj^^(‘^)i ^or weak scattering by a power-law 
phase screen. These spectra were computed numerically from equations 
A38 and 2.4. The phase alternates between values of 0 and n, with 
abrupt shifts occurring at frequencies « = (m) ' for Integer values 

of m. Each shift In phase Is accompanied by a null In the magnitude 
of w *(<*>); the depth of the nulls varies as a result of discrete 
sampling. For these calculations, e/f = 10, |v*j = 7.07 km s , 

Y = 90°, A = 0°, p = 3.5, (XO)^^^ = 5 km, and w = w' = 2n. 

W Va 
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APPENDIX D 


THE EFFECT OF THE QUIESCENT ATMOSPHERE DURING PLANETARY 

RADIO OCCULTATIONS ' 

The effect of the quiescent atmosphere during planetary radio 
occultatlons has been discussed at length in the literature (for 
example, see Haugstad [1982]); this appendix summarizes the results and 
defines parameters. Additional references on this subject include 
Fjeldbo and Eshleman [1969], Young [1976], Haugstad [1978a, ll'Bc, 

1982], Hubbard et al. [1978], Eshleman et a1. [1979a], and Woo et el. 
[1980b]. 

Figure D1 represents the idealized occultation geometry and shows 
the coordinate system. A monochromatic plane wave propagating in the x- 
direction impinges on a planetary atmosphere which is free of smell- 
scale irregularities such as turbulence, waves, or layers. The 
refractive index In the atmosphere is given by 1 + n^(r). Afte.’ 
refraction by angle a, where <* << 1, the geometrical optics ray which 
Intercepts the spacecraft has fractional Intensity where the two 

factors account for focusing due to limb curvature and deforusing due to 
the vertical gradient in n^(r), respectively. For a spherically 
symmetrical. Isothermal atmosphere with scale height H, 

= and = . (Dl) 

Here, D is the distance from spacecraft to planet, and r Is the altitude 
of the raypath at periapsis. 
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Effect of the Quiescent Atmosphere 


As a second consequence of differential refraction, the transverse 
raypath velocity in the phase screen is given by 


■* 



Ij) y u 

y y' 




(D2) 


where v^ and v^ are components of spacecraft velocity that are 

* A 

perpendicular to the direction of propagation, and y and z are the 

horizontal and vertical coordinate unit vectors at ray periapsis. 

Finally, the quiescent atmosphere distorts the Fresnel zone at 

1/2 

wavelength A, to horizontal and vertical dimensions (<PyAD) ' 
and respectively. 



S/C 


Figure 01 Idealized occultation geometry and coordinate system. A 
radio signal transmitted from earth (located at x “ undergoes 
refraction by angle ut in the f.cculting planet's atmosphere and is 
received subsequently by the spacecraft. A distance 0 separates the 
spacecraft and planet; the parameter r gives the altitude of the raypath 
at periapsis. As viewed from earth, the y- and z-directions are tangent 
and perpendicular, respectively, to the visible limb of the planet. 
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APPENDIX E 


THE VALIDITY OF THE MEASUREMENT UNCERTAINTIES 

The results reported in Chapters 4 and 5 consist of solutions for 
four parameters that characterize the refractive irregularities In the 
Ionospheres of Jupiter and Saturn. From Chapter 6 It Is clear that 
realistic uncertainties are essential for Interpreting these parameter 
measurements. In this appendix, the uncertainty In y, the inferred 
magnetic field orientation (see Chapter 6), is shown to be realistic 
with respect to statistical errors and not simply formal. As the 
uncertainties In all four parameters were calculated by the same 
procedure (Section 3.4), this suggests, by implication, that the 
uncertainties in the remaining parameters are also reliable. 

The conclusions In this appendix do not apply to systematic errors, 
which are more difficult to assess. One such error in y arises from 
uncertainty in the spacecraft trajectory. This effect was considered, 
but was found to be Insignificant for the occultations of Voyager by 
Jupiter and Saturn. As mentioned previously, a bias in r to the 
spatially varying field geometry in Jupiter's ionosphere could be 
comparable to the statistical errors; however, this source of error was 
taken into account in evaluating the Jovian magnetic field models 
(Section 6,3). The analagous effect at Saturn is believed to be 
negligible relative to the statistical uncertainties (Section 6,4). 
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Validity of the Measurement Uncertainties 


The 13-cm amplitude spectrum measured 1n region 2 during V2 
emersion at Saturn (see Figures 5.4 and 5.5) Is used as an example In 
the first part of this discussion. In order to demonstrate the 
sensitivity of the amplitude specti o changes In y. three theoretical 
spectra are shown In Figure El. One curve (y = 46°) corresponds to the 
least-squares solution for parameters derived from this measured 
spectrum (see Table 5.3); this same theoretical spectrum also appeared 
in the panel at the lower left of Figure 5.5. The other two curves 
Illustrate the effect of changing y by ±3o, or 9® (see Table 5.3), from 
the best fit. Note that a shift of 3a in either direction significantly 
alters the alignment of the theoretical spectrum, a phenomenon which was 
explained In Appendix C. 
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Validity of the Keasuremijnt Uncertainties 



Figure El Theoretical 13-cm amplitude spectra for weak scattering 
by a power-law phase screen. The spectra were computed numerically as 
described in Section 3.4. One curve (y = 46®) corresponds to the least- 
squares solution for parameters derived from the 13-cm amplitude 
spectrum in region 2 during V2 emersion at Saturn (see Figures 5.4 and 
5.5); this same theoretical spectrum also appeared in the panel at the 
lower left of Figure 5.5. The other two curves illustrate the effect of 
changing y by ±3a, or 9“ Jsee Table 5.3), from the best fit. Note that 
a shift of 3o in either direction significantly alters the alignment of 
the theoretical spectrum, as explained in Appendix C, For these 
calculations, e/f = 7, |v*| = 20.6 km s p = 3.6, and 
= 5.21 km, where X is the radio wavelength. 
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Validity of the Measurement Uncertainties 


In Figure E2, three theoretical curves are compared to the measured 
amplitude spectrum. For these fits, the angle y was assigned a fixed 
value {37®, 46°, or 55°, as in Figure El) and the three parameters p, 
e/f, and Q (see Section 2.3,1 for definitions) were adjusted in order to 
minimize the residual, as defined in equation 3.5. The residual attains 
its overall minimum value when y = 46°, as reported in Chapter 5 (see 
Figure 5.5 and Table 5.3). A change in y by ±3a away from the least- 
squares solution visibly reduces the quality of the fit in the frequency 
interval from 1 to 4 Hz. (The solutions for p, e/f, and Q for each 
choice of y remained within la of the values listed in Table 5,3.) 


162 


Validity of the Measurement Uncertainties 



FREQUENCY, Hz 


Figure E2 Comparison of three theoretical curves to the 13-cm 
amplitude spectrum measured in region 2 during V2 emersion at Saturn. 

For these fits, the angle y was assigned a fixed value (37°, 46°, or 
55°) and the other three parameters in the theory (p, e/f, and Q) were 
adjusted in order to minimize the residual. The overall minimum value 
of the residual occurs when y = 46°, as reported in Chapter 5. A change 
in y by ±3a away from the least-squares solution visibly reduces the 
quality of the fit in the frequency interval from 1 to 4 Hz. 
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Validity of the Measurement Uncertainties 


The residual of the theoretical fit to the measured amplitude 
spectrum was calculated as a function of y (by varying p, e/f, and Q) 
for values intermediate to those shown in Figure E2. The results appear 
as curve 1 in Figure E3; the angle y has been converted to r (see Figure 
6.1) for comparison with the predictions of magnetic field models. The 
residuals of curve 1 at r = 217®, 226®, and 235® correspond to the three 
fits in Figure E2; the well-defined minimum at r = 226® indicates the 
le«ist-squares solution for this parameter. When r is shifted by 
±3<J away from the minimum, the residual increases by more than 70%. The 
results in Figures E2 and E3 demonstrate that the occultation 
measurements are sensitive to changes of 3o in the magnetic field 
direction, and confirm the uncertainty in y- 

Figure E3 also shows plots of the residual as a function of r for 
the 13-cm phase spectrum from region 2 (curve 2) and the 3.6-cm 
amplitude spectrum from region 3 (curve 3) of V2 emersion at Satur>i (see 
Figures 5.4 and 5.5). The corresponding curves derived from the two 
spectra measured during VI immersion at Saturn (see Figures 5.2 and 5.3) 
appear in Figure E4. The similarities among all five curves suggest 
that the amplitude spectrum discussed above is typical of the Voyager 
occultation measurements. 

For comparison. Figures E3 and E4 also indicate the values 
of r predicted by the OD and Z3 magnetic field models (see Section 
6.4). The minimum of the residual in each of the five curves is 
intermediate tu the model predictions, as discussed in Chapter 6. 
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Validity of the Measurement Uncertainties 



Figure E3 Normalized residual as a function of r. Curve 1 was 
derived from the 13-cm amplitude spectrum from region 2 of V2 emersion 
at Saturn (see Figures 5.4 and 5,5). To calculate one point in the 
curve, r (or equivalently y) was assigned a fixed value and the other 
theoretical parameters (p, e/f, and Q) were adjusted In order to 
minimize the residual. The values at r = 217°, 226°, and 235° 
correspond to the three fits In Figure E2; the minimum at 226° indicates 
the least-squares solution for this parameter. When r is shifted by 
±3a away from the minimum, the residual increases by more than 70%. The 
curves corresponding to the 13-cm phase spectrum from region 2 (curve 2) 
and the 3.6-cm amplitude spectrum from region 3 (curve 3) of V2 emersion 
at Saturn are also shown, along with the angles predicted by the 00 and 
Z3 magnetic field models. 
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Validity of the Measurement Uncertairties 



164 169 174 

r, degrees 


Figure E4 Normalized residual as a function of r for the 13-cm 
amplitude spectra measured in regions 1 and 2 of VI Immersion at 
Saturn. See Figure E3 for additional comments. 
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Validity of the Measurement Uncertainties 


An empirical test of the uncertainties comes from the measurements 
of Y In Tables 4, 2-4. 4 and 5. 2-5. 3. These solutions for y were obtained 
either from non-overlapping time segments of data, from phase and 
amplitude spectra measured over the same interval, or from measurements 
at different radio wavelengths, so that all the measurements at any one 
location are based on different data sets. As the magnetic field 
orientation at Jupiter or Saturn Is not expected to vary appreciably 
with tltltude over the Ionospheric region of the observations, these 
measurements should be mutually consistent at any one location. Of the 
19 measureinents of y at Jupiter and Saturn, 16 are within lo and all are 
within 2.5a of the average value of y at their respective locations. 

This result supports the accuracy of the uncertainties. 
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